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ABSTRACT 



Many Navy applications of composites , , including ships' superstructures, 
submarine air flasks and missile rocket motor casings, require high 
strength and reliable materials. Composite strength reliability is dictated by 
individual fiber breaks at low loads (lower tail) and the accumulation of the 
fiber failure sites. 

This study examined the effects of applying a preload to a graphite/epoxy 
composite tow prior to complete polymerization of the matrix. The 
objective was to break the (inevitable) weak fibers and minimize the effects 
of the associated stress concentrations, subsequently limiting the clustering of 
fiber failures. By eliminating the lower tail, the shape of the Weibull 
distribution is reduced, thereby enhancing composite reliability. 
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I. INTRODUCTION 



In the 1960s, Composite Materials were referred to as "Space Age 
Materials". Today, composite structures are common through out our 
society. The rapid growth of the composite materials industry has been 
phenomenal, perhaps only exceeded by the advancements in computer 
technology. The concept of composite materials is simple, in that two or 
more structurally compatible materials are combined to produce an end 
product with attributes greater than the sum of the constitutes. But as new 
applications for composite materials are developed so follows an increase in 
unanswered questions with regards to reliability, especially associated with 
vital national defense systems, public transportation systems and aerospace 
and hydrospace exploration. As engineers it is important that we do not let 
the state of the art exceed the state of the science. We must have answers to 
those difficult questions concerning composite material applicability, 
maintainability and reliability. 

Composite reliability is the focus of attention at the Advanced 
Composites Laboratory at the Naval Postgraduate School, Monterey 
California. Here questions are being investigated regarding composite 
reliability as a function of service life, strength and test methodology. The 
primary goal of composite material research is to be able to have a better 
understanding of how composite materials behave over an extended service 
life and to further investigate the mechanics of composite strengthening. A 
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secondary function is the investigation into improved testing methods so that 
research techniques can be translated into production applications. 

The primary objective of this study was to investigate composite 
reliability as a function of the strengthening mechanisms. This was 
accomplished by studying the affects of preloading a graphite bundle prior to 
complete curing of the epoxy matrix. After the matrix had cured, the 
graphite strand ultimate tensile strength was determined and the results were 
compared to a benchmark value based on previous testing. 

In addition, individual fiber testing was conducted so that a 
mathematical model may be developed in order to determine the optimum 
bundle preload level. Fiber testing was also conducted in order to 
demonstrate that the Integrated Fiber Testing system developed at the Naval 
Postgraduate School (NPSIT) could be use as a viable and accurate test 
method. This was accomplished by conducting fiber testing using an 
INSTRON Model 4200 materials tester using current American Society of 
Testing and Materials (ASTM) procedures and comparing the results 
obtained using the NPSIT system. Graphite fibers having different strength 
properties were tested in order to show that the testing systems used could 
accurately ascertain any small differences in composite properties. The 
focus of this study was on fiber diameter and fiber tensile strength. 
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II. BACKGROUND 



Load bearing graphite composite materials consist of high strength / high 
modulus graphite impregnated with a ductile reinforcing matrix. The 
performance characteristics of the composite material are controlled by 
three factors: 

- Strength and modulus of the graphite fiber 

- Strength and modulus of the epoxy matrix 

- Effectiveness of the bond between the matrix and fiber, with 
regards to the mechanics of load transfer at the interface. 

The overall tensile failure process of this simple system is extremely 
complex, however for the purpose of this study this can be simplified by 
assuming a two dimensional model and that failure initiation occurs within 
the fiber only. This failure model was first presented by B.W. Rosen (Ref.l 
and 2). It consisted of parallel fibers in a homogeneous matrix, loaded in 
tension, assuming a uniform strain supported primarily by the fibers. 

The most important concepts developed in this classical study were the 
fact that the fibers were considered to have statistically distributed flaws or 
imperfections. Theses flaws resulted in fiber failure at various stress levels, 
therefore Rosen considered each fiber as a series of links and that the 
individual fiber failure process was a function of the weakest link principle. 
Finally, that the statistical strength distribution of theses flaws could be 
approximated by a Weibull distribution. 
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A. FIBER FAILURE PROCESS 



Initially when a fiber bundle .without epoxy, is loaded in tension, the 
individual fiber failures occur randomly at the site of a flaw at a load much 
less than failure load ( x\ < Xf ). When this occurs the stress (01) in the fiber 

becomes zero, and the stress in the surviving fibers increase to a new stress 
level equal to some value K * cq. The existence of some individual fiber 

failures may not be catastrophic in that the increase in stress is shared equally 
by the remaining fibers as long as K * cq < Of. Then applied tension x can be 

increased to a new value X2 > x\ then this process is repeated until * cq > 
af. Excluding any effects due to twist or friction, catostophic bundle failure 

occurs after each individual fiber breaks once. 

When the fiber bundle is impregnated with an binder matrix the 
mechanical behavior of the composite strand is altered. The addition of the 
matrix creates a unique load transferring system. This provides for 
transverse load sharing between each individual fiber and longitudinal load 
sharing to another segment of the fiber. The matrix also provides the 
mechanism for localizing the effects of microcracks within the composite 
material. These two functions combined provide for the distinct mechanical 
performance of fiber reinforce composite materials. 

The composite failure process based on Rosen's failure model, with 
matrix binder, predicts different results from dry bundles. When the initial 
load Xi is applied initial fiber failures still occur at the inevitable flaw sites, 

however the effects are confined to a localized region. Stress o decreases to 
zero, and the shear x in the matrix reaches a maximum value. The axial load 
is transmitted by this shear to only the neighboring fiber segments, (see 
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Figure 1). The neighboring fibers now show a load concentration due to the 
fiber break. At the broken fiber tip there exists a length of fiber which is 
ineffective in carrying the applied load. This is known as the ineffective 
length 5. The increase in shared stress on the immediate neighbors is 
substantially greater then that shared by the fibers that are a greater distance 
from the failure site. 



t 

Load 
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Figure 1. Fiber Bundle Model 
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The determination of the ineffective length 5 is difficult, and requires an 
understanding of the shear stress distribution along the interface. The 
ineffective length is estimated by Rosen [Ref. 1] to be: 

5 = 1/2 {( Vf 1/2 - 1) E f /G m } 1/2 cosh-1 {[1 + (1 + <p)2] / 2[1 - q>] } d f 
Vf - Volume fraction of fiber 
Ef - Modulus of fiber 
Gf - Shear modulus of matrix 

(p - Fiber efficiency or the fraction of the undistributed stress 

value below which the fiber is considered to be ineffective 
df - Fiber diameter 

B . COMPOSITE FAILURE 

In general there are three basic modes of composite failure. The first is 
debonding which is failure at the fiber/matrix interface which is caused by a 
high interface shear stress (Figure 2a). This failure may start with a single 
broken fiber and propagates within the interface, along the fiber length, 
which in turn drastically increases the ineffective length 5. Thus increasing 
the likelyhood of the crack propagating to another flaw site. The study of 

this type of failure identifies the need for improving the interface properties 
such as adhesion and toughness properties of the epoxy (Note: xj « Xf). The 

second type of composite failure is crack propagation (Figure 2b). This 
occurs when an initial crack is developed at the flaw site and propagates 
transversely across the composite. Crack propagation is controlled by the 
fracture toughness of the matrix and fiber and thus minimizing the effects on 
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local load sharing. This failure mode may be associated with some slight 
debonding. (Note: Xj « Xf). 




Figure 2a. Figure 2b. Figure 2c. 



Figure 2. Composite Failure Model 

Finally, if debonding and crack propagation are controlled then 
composite failure occurs when as load is increased, the number of statistically 
distributed fiber failures accumulates in the vicinity of one cross section such 

that when the load shared in the surviving fibers exceeds the ultimate load of 
the system ( xj > Xf ). This mode of failure is usually a combination of 

debonding and crack propagation, occurring at a much higher load level, 
(Figure 2c). 
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C. LOAD SHARING 



The most pronounced effects of adding a matrix binder to the fiber 
bundle is that of load sharing. This is best discussed using a microbundle 
model developed by S.L. Phoenix. [Ref. 3]. This model, (Figure 3) is made 
up of n fibers, each consisting of m links having a length of 5, the ineffective 
length. In each microbundle ( = n * 5) there exists a set of statistically 

distributed flaws. When the load x is applied the fibers will sequentially fail, 
in accordance with a Weibull distribution function. 

F(x) = 1 - exp [-(x/(3 x ) a ] 

Where a is the shape parameter and p is the scale parameter. 

Because the individual fiber consists of m links , the weakest link 
theorem applies. Therefore, the value of p x is a function of gage length and 

must be corrected for GL. 

P2x = Plx (GL 1 /GL 2 ) 1/a 

When an initial fiber failure xj occurs in a weakest link the load on each 
immediate neighbor increases to a value X=Ki*x, where Kj is the load 
concentration factor and Kj> 1, (Figure 3). If another weak link exists within 

the region of concentrated load, then a second fiber failure will prematurely 

occur. This in turn creates an additional load sharing burden on the 
neighboring fibers X = K 2 *x, where K2 >Kj> 1, (Figure 3). Therefore, 

randomly created load concentrations may prematurely break a larger 
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number of fibers at a load xj than would otherwise normally occur, 

subsequently weakening the composite. J.M. Hedgepeth [Ref. 4] provided 
quantitative values to the load concentration factors, where K\ = 1.33, K 2 = 

1.60, K3 = 1.83 and K4 = 2.03. Meaning, for example, that the load in the 
fibers adjacent to three consecutive fiber breaks equals 1.83 * applied load . 



t 

LOAD 




If the random fiber breaks occuring at Xj were dispersed throughout the 

composite the effects of load concentrations K would be minimal, however 
this can not be controlled. Therefore an attempt must be made to control the 
effects of the spatial clustering of fiber breaks and load concentrations in a 
local region. 
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A.S. Tetelman [Ref. 5] and C. Zweben [Ref. 6] provide excellent 
summaries of the composite failure process and the effects of load 
concentrations. Tetelman also discussed a few concepts for minimizing load 
concentrations and controlling crack propagation. These include the use of 
discontinuous fibers, using fibers or matrix consisting of various modulus of 
elasticity and the use of prestressing at elevated temperatures In fact 
prestressing composite bundles at elevated temperatures was proven to be a 
viable solution during experimental work conducted by E.M. Wu. [Ref. 7]. 
G.J. Mills [Ref. 8] has also studied the effects of prestressing Boron/Epoxy 
prepregs by rolling. 

D. COMPOSITE RELIABILITY 

The reliability of a composite structure depends on a definite 
characterization of the weak strength distribution (the weak lower tail). 
Definitive evaluation of the lower tail requires a large sample size and 
consistent material production system. 

Such a data base for composite material performance is usually 
developed from either limited fiber/bundle testing or complete system testing 
of a small sample size. A small data base, when combined with the existence 
of the statistically distributed flaws or imperfections in the graphite fiber 
due to manufacturing, processing or handling practices, makes the behavior 
of a given composite structure difficult to predict. When a multimillion 
dollar rocket motor casing or aircraft component consists of an enormous 
number of fiber links, performance reliability becomes dependent upon a 
statistical study of the strength characteristics. Because of the existence of 
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the fiber flaws, the strength values of the Weibull distribution are well 
dispersed within the lower tail. This imposes great restrictions on the 
allowable design limits for the composite structure. Ideally, there would 
exist nearly no dispersion in the test data (coefficient of variance approaches 
zero) and then the data would be completely reliable for a given strength 
level. 

This study will focus on improving composite strength reliability by the 
use of preloading the fiber bundle prior to the curing of the epoxy resin in 
order to remove the weak fiber sites (lower tail of the Weibull distribution) 
and decrease the dispersion of the composite strength. 
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HI. EXPERIMENTATION 



The experimental work was concentrated in two areas. First, 
determining the diameter and strength of the individual fiber and second, 
determining the strength characteristics of a graphite tow subsequent to 
preloading. 

The material tested was from a Hercules Magnamite high strength 
graphite, type AS-4 spools 008 and 019. The graphite bundle consisted of 
3000 fibers with a nominal diameter of 7 micrometers and a denier of 
.005746 grams/inch. 

A. FIBER TESTING 

Individual fiber testing was a vital part of this study which provided the 
essential data required to determine the appropriate preload. First the fiber 
diameter was determined then the fiber failure strength. This information 
was used to develop a Weibull probability plot to statistically determine the 
number of fiber failures as a function of applied load. 

To study the feasibility of using an alternative testing method, two 
testing methods were used during this study. The first method was the use of 
an Instron Material Tester (INSTRON), (Figure 4), with associated ASTM 
procedures to determine fiber failure load. Second, using a Integrated 
Testing system under development at the Naval Postgraduate School, 
Monterey, California. (NPSIT), (Figure 5). 
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Figure 4. INSTRON Testing System 




Figure 5. NPSIT Testing System 
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The primary advantage to the NPSIT system is that the fiber diameter as 
well as the failure load can be readily determined on the same test stand. This 
minimizes the handling of the test samples and bias the statistical parameters. 

The correlation between fiber test length and failure load has been well 
established in past studies and is not addressed here. A gage length of 5 centi- 
meters was used for fiber testing. Strain rate or crosshead speed also has an 
effect on test results, however this effect was not studied. A constant speed of 
10 percent of gage length per minute (in accordance with ASTM procedure) 
was employed. 

Regardless of the testing method used, sample preparation (Appendix 
A.) and testing parameters were identical. 

1. Fiber Diameter 

Current methods for determining fiber diameter, as outlined in 
ASTM procedures, call for the use of the Scanning Electron Microscope 
(SEM). This is a time consuming and some what arduous task when testing 
large numbers of fibers. 

Recent work at the Naval Postgraduate School, Monterey, has provided two 
additional methods which are more conducive to this study. The principles 
and procedures as developed by T. A. Bennett [Ref. 9] are based on the 
physical optics, that if a illuminated beam of light is directed at the fiber that 
a diffracted light pattern is produced which is perpendicular to the fiber 
( Figure 6). 
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Figure 6. Diffraction Pattern 



The first method Bennett looked at was the use of a 
photoconductive cell to measure the distance between nodes (null intensity in 
the diffracted light pattern). The second method consists of the use of a 
microneye interfaced with an Apple Plus II microcomputer. Bennett [Ref. 9] 
and M. Storch [Ref. 10] then developed the data acquisition techniques 
necessary to accurately determine fiber diameter. The required test 
configuration is essentially the same for both methods, however, since only 
the photocells was used in this study , only those associated procedures will be 
further discussed. 

The test equipment ( Figure 5) consists of ; a low power (.052 mW) 
Helium-Neon laser, which provides the required light source of known wave 
length ( lambda = 632.8EE-09 meters). A focusing lens, a spatial filter and a 
collimating lens are enclosed as one unit and mounted on the end of the laser. 
This is necessary in order to produce the parallel light pattern. The fiber is 
mounted on the test stand which by use of micrometers allows for 
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longitudinal and transverse adjustments, rotation of the sample and 
adaptability to the load testing device. At the end of the tracks are two 
photocells mounted to adjustable pedestals which allow for transverse 
adjustments. The photocells are connected to two Fluke 8840A multimeters. 

2. INSTRON Load Test 

Individual fiber load testing was performed using the Instron Load 
Tester Model 4206 (INSTRON) inconjunction with compatible data 
acquisition systems. The data acquisition system to be used during this study 
is the Hewlett-Packard 3497A Data Acquisition/Control Unit and the 
Hewlett-Packard HP-85 for computer programming. Appendix B provides 
a detailed listing of the Interactive Data Acquisition Software (IDAS). IDAS 
was design to minimize operator effort and ensure accurate and 
reproducible test results. All IDAS procedures were written in accordance 
with American Society of Testing and Materials (ASTM) Standard 
Specifications. 

After the required program is loaded enter RUN and then follow 
the outlined steps in Appendix B. The programs use the function keys to 
provide for flexible data acquisition. Options are normally provided for a 
CRT/printer output or with the graphics plotter. 

Fiber diameter was measured prior to this test with the use of the 
NPSIT test system. 

The IDAS software includes the following programs. 

LDCALB - Load Cell Calibration: Characterization of the 

strength of a fiber required absolute calibration of the load transducer. This 
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was accomplished by calibrating the digital output of the load transducer by 
standard weights at six levels. This provided absolute calibration all 
analogue and the analogue to digit interface. 

MCTST - Characterization of the stiffness of a fiber requied the 
calibration of the system compliance. The total compliance measured (C t ) is 

the sum of the fiber compliance (Cf) plus the system compliance (C s ), all 
three are a function of applied load level P. 

Q(p) = ckp) + c s ( P ) 

The system compliance C s is comprised of the compliance of the 
testing machine, of the cardboard tab and the adhesive which bonds the fiber 

to the tab. If the system compliance is assumed to be constant, then a single 
calibration constant ( effective gauge length of zero) is determined as 
recommended by ASTM procedure D-3379-75(82), (Ref. 11). as discussed 
in ASTM procedure 3379-75(82), the effective gauge length can be 
extrapolated from measurement of the total compliance using three different 
gauge lengths. The short coming of this procedure is that it assumes that 
system compliance is independent of load level P. The latter was observed to 
be a gross assumption. Compliance calibration was accomplished using an 
actual zero gage length sample. In this case, the fiber compliance Cf term 

vanishes and the total compliance measured is the system compliance. 
Furthermore, the compliance, over the entire load range of the strength of 
the fiber measured was recorded and the data fitted by a second degree 
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polynomial. This second order polynomial was then used to calibrate the 
tensile strength test data, thereby compensating for the non-linear depencey 
on the load level. 

QFIT - Quadratic Curve Fit: This program reads the saved 

Machine Compliance data and provides the three coefficients of a quadratic 
curve. These values E, F, G were then recorded for further use in the fiber 
test program to mathematically subtract from the total displacement the 
displacement associated with the system compliance. 

FTST - Fiber Test: Is the primary program in this series. 
INPUTS, allows for the input of all key parameters necessary to determine 
the actual fiber ultimate strength. The test is run for 30 seconds at a speed of 
.lx gauge length (5 millimeters/minute on a 50 millimeter sample). The A, 
B, C, E and F inputs as discussed above are required to perform the necessary 
mathematical calculations. AQUIR, provides for the actual proof test and 
data acquisition. Once the sample is ready the program will adjust for the 
initial tension then start the timer, providing for 300 data readings during the 
30 second, run. The results can be displayed on the CRT/printer using GRAF 
or on the HP Plotter using PLOT. SAV-DAT and RD-DAT provide for the 
storage and reading of test data. CAL, mathematically subtracts system 
compliance from the total displacement and provides the unbaised results. 
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LDTST - Load Test: Provides for a quick check of the load cell 
calibration using the original A, B and excitation voltage from the LDCALB 
results. 



LDSTOR - Load/Store: is a data management tool for transfering 

files. 



3. NPSIT Load Testing System 

The Naval Postgraduate School Integrated Load testing system was 
developed as an alternative testing method which provides for the 
determination of fiber diameter as well as failure load with one test system. 

The data acquisition and software is essentially the same as that used 
with the INSTRON test system. The only exception being that a short 
program was added to FTST so that fiber diameter and failure load are 
recorded on the same file for easier access. 

By following the procedures in Appendix A and Appendix B, the 
NPSIT system can be quickly initialized and ready for fiber testing. After 
the FTST program is loaded, mount the fiber sample in the clamps. First 
measure the fiber diameter, as discussed and then enter AQUIR. The distinct 
difference between this test method and the INSTRON test method is that the 
NPSIT drive motor is manually controlled by the operator upon prompt 
from the HP-85 program, where as the INSTRON tester is software 
controlled. 

All fiber data is recorded, displayed and stored in the same manner 
as with the INSTRON test. 
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B. TOW TESTING 



The objective of graphite tow testing was to determine the the effects on 
the strength characteristics as a result of applying a preload. Experimentally, 
AS-4-019 graphite bundles, from the same spool used in fiber testing, were 
tested. 

1. Initial Set-up 

The samples were prepared by paying the graphite out on a table 
and applying a 2 kilogram load. Then copper tabs (1 inch x 1.5 inch) were 
fasten to the bundle at a gage lengths of 10 inches using an epoxy mix as the 
adhesive. The 2 kilogram load removes the slack and aligns the individual 
fibers within the bundle. The copper tabs provide for an area for the bundle 
to be gripped in the Instron tester. 

The Instron Universal Testing Instrument Model 4206 was used to 
provide the tensile load to the test samples. This system was used in 
conjunction with the 4200 Series Expanded Control Console with data 
acquisition being provided by the use of an IBM PC-AT. The Instron/IBM 
data acquisition on system provides for the real time graphical display of load 
versus displacement and the creation of a data base for additional 
computations and graphical outputs. The initial operator and report files are 
established prior to the commencement of testing to generate the desired 
reports and graphs. 

2. Preloading Bundles 

Prior to applying a preload, the Graphite bundle was impregnated 
with the epoxy matrix. The matrix consisted of a mixture of Dow DER-332 
resin (55 percent by weight) and Texaco Jeffamine T403 (45 percent by 
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The final resin content of the strand was determined by weighing a 
few random samples and comparing to the weight of the dry bundle. 

% R = W T - W P 
W T 

Where Wj is the total weight of the strand and Wp is the weight of the 

bundle (denier x gage length). A resin content of 55 to 60 percent is desired. 
3. Load Testing 

The failure load of the final product (the matrix impregenated 
composite strand) was determined using the Instron/IBM system described 
earlier. After the strand is mounted in the test grips the desired test files are 
inputted and the machine crosshead speed is set a 10 centimeters per minute. 
The test is allowed to continued until complete tensile failure of the 
composite strand occurs. The final output provides for the maximum load in 

kilograms and displacement at failure in millimeters. A graphical report 
may also be generated. 
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VI. RESULTS 



A. FIBER TESTING 
1. Fiber Diameter 

A review of the statistical summary of the fiber diameter results 
(Table I), provides a clear indication that the test method used for 
determining fiber diameter produced consistent and accurate results. The 
normal probability plot, Figure 8, shows how the AS-4-008 series and the 
AS -4-0 19 series graphite diameters compared to the merged results. 

A total of 79 samples were measured from the AS-4-008 graphite 
spool, with a mean diameter of 7.233 microns versus 82 fiber samples from 
the AS-4-019 spool with a mean diameter of 7.245 microns. The large 
variability of fiber diameter (as measured by the magnitude of the standard 
deviation) verifies that the fiber diameter is not consistent over the entire 
length. 

Exploratory effort was expended in the adjustment of the test stand, 
optimizing the distance between the fiber and the projected plane (the 
photocell).. This was required to compensate for the narrow width of the 
photocell window. It was necessary to adjust the test stand so the width of the 
node was nearly as wide as the photocell. If the node is too wide, there exists 
several locations of maximum resistance readings. If the the node is too 
narrow, too much light saturates the photocell, thus making it impossible to 
locate the center of the node. 
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008 Series N Shape Scale Mean Std Dev Coef. 

Load (gmf) OC j3 Var. 



TABLE I: STATISTICAL SUMMARY OF SINGLE FIBER TESTING 
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Figure 8. Normal Distribution Plot of Fiber Diameter 
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2. Fiber Failure Load 



Interpretation of fiber failure data is focused on two areas 

— Comparison of test results from INSTRON and NPSIT testing to 
identify any differences. 

— Comparison of the statistical strength for samples from two different 
spools. 

In Table II statistical parameters from Table I are employed for 
comparison of the INSTRON and NPSIT test methods. (see Figures 9 and 10). 

TABLE n. FIBER LOAD RESULTS 
Mean(gmf) Std Dev 





008 


019 


008 


019 


INSTRON 


14.26 


13.91 


4.20 


4.82 


NPSIT 


14.56 


14.06 


2.78 


3.20 


Difference 


+2% 


+1% 


-15% 


-15% 



Comparing the first two columns of TABLE II, it is apparent that the mean 
strengths of the samples measured with the NPSIT are, slightly but 
consistently, higher for both spools tested. Comparing the last two columns, 
the variabilities of the measured samples in the NPSIT are, significantly and 
consistently, lower for both spools. These observation substantiate the 
hypotheses that the NPSIT, which integrates the fiber diameter measurement 
with the strength measurement, enables the elimination of several handling 
and storage steps, thereby minimizing the associated accidental damage. 
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Figure 9. Weibull Distribution Plot of Fiber Load 
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Figure 10. Weibull Distribution Plot of Fiber Load 
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Minimization of this damage results in a shift of the lower tail, of the Weibull 

plot, to the left (evident in both Figure 9 and Figure 10 for samples labeled 
*) and a significant reduction of variability and the standard deviation. 

Minimization of the handling damage cannot increase the instrinsic strength, 
therefore has no effect on the upper tail (again evident in Figures 9 and 10). 
The shift of the lower tail changes the geometric centroid, of the probability 
distribution function (pdf), therefore causing a slight increase of the mean. 
The consistency between the quantitative measurement with the qualitative 
expectations support the conclusion that the integrated NPSIT test is a worth 
while improvement. 

For comparison of the statistical strength of the two spools, the 
strength is modeled by the two parameter Weibull distribution, Table III. 
The parameters, as based on Maximum Likelihood Estimation, of the 
respective data sets are summarized: 

TABLE HI. FIBER LOAD WEIBULL PARAMETERS 







Alpha 


a 




Beta p 






NPSIT 


INSTRON 


MERGE 


NPSIT 


INSTRON 


MERGE 


008 


5.75 


3.62 


4.28 


15.72 


15.80 


15.79 


019 


4.86 


3.32 


3.94 


15.33 


15.52 


15.45 


Diff 


-18% 


-9% 


-8.5% 


-2.5% 


-1.8% 


-2.2% 



Again, note the significant difference between the a parameter 
(which is inversely associated with variability) for the two spools. This is 
particularly evident when handling has been reduced On the other hand, 
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consistent with the above discussion, the difference between the p parameter 
(associated with mean strength) is slight. Physically, composite samples 
fabricated with the 019 spool which has a smaller a (or larger variability) 
will have a larger number of fiber breaks than those fabricated with the 008 
spool. In accordance with the local load sharing model, the larger number of 
fiber breaks will increase the stochastic process of clustering of failure sites 
hence, leading to a higher composite scatter. Therefore, preloading spool 
019 and observed accompanied change of composite variability will offer 
both a verification of the load sharing model and explore the practicality of 
improving composite reliability. 

Complete tabulated results for all fiber diameter and tensile 
strength measurements are found in Appendix D. 

B. TOW BUNDLE PRELOAD 

The objective of preloading the tow bundle, prior to matrix binder 
impregnation was to decrease the variability of the composite strength. It 
was established in [Ref.l] and [Ref. 2] that composite failure occurred 
sequentially when weak fiber failure began at low load levels. As applied 
load is increased, the number of broken sites increased and formed spatially 
clustered failure sites leading to stress concentrations, the most severe of 
which causes catastrophic failure. This failure process can be observed by 
comparing the load deformation curve of a matrix-free tow bundle (curve B 
in Figure 11) and a matrix impregnated bundle, the composite (curve C). 

For the tow bundle curve B, the initial slope of the load deformation 
curve is the sum of all the filaments in the bundle. As the deformation 
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Figure 11. Composite Deformation Curve 
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Percent (Z) 



displacement is increased, the weakest fibers break and the load carrying 
capability of the entire bundle is lost, thus resulting in a decrease in slope 
such that the secant modulus is equal to the sum of modulus of the surviving 
fibers. This continues until only a few very strong fiber survive and the 
secant of the load deformation curve of the bundle approaches horizontal. 

For the composite curve C, the initial slope of the load deformation 
curve is the same as that of the bundle curve B, since both have the same 
number of filaments. Within the composite, weak fibers start to fail at 0.3 
percent strain as in the bundle. However, in the presence of the matrix 
binder, the load carried by the broken fiber is transferred to the neighboring 
fibers via the matrix and, instead of the entire fiber, only a portion equal to 
the ineffective length is lost. The ineffective length is approximately equal to 
10 fiber diameters or in the case of graphite, in the order of lOOum. Such a 
minute loss can not be detected from the resolution of the load deformation 
curve. As result the load-deformation curve remains apparently linear up to 
the point of catastrophic rupture even though numerous filament failures 
sites accumulated internally. 

Therefore, it is desirable to develop a preload procedure to break the 
inevitable weak fibers, without creating stress concentrations which lead to 
premature undue breakage of the neighboring fibers. This is possible if the 
preload can be performed before the matrix can effectively transfer the load. 
Subsequently when the matrix becomes effective, the load around the broken 
site can be transferred to the neighbor without causing spatial clustering of 
broken fiber sites and subsequently reducing the variability of the ultimate 
strength. 
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The most obvious implementation is to perform the preload on a dry 
bundle prior to the introduction of the matrix binder. Exploration proved 
that this approach was not practical. In the absence of matrix binder, a 
broken filament unleashed the stored elastic energy dynamically leading to 
fraying and lost of alignment of the filaments, (see Figure 12). In tact this 
phenomenon occurred when preload was applied to a bundle wetted m a 
matrix but prior to polymerization. In such a case the surface tension of the 
uncured polymeric matrix is insufficient to contain the dynamic energy 
release of the breaking filaments. Therefore, the impregnated bundle was 
allowed to air cure pnor to preloading to ensure sufficient surface tens.on to 

prevent fraying. 
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Based on the results from Appendix C, the preload level at a 
deformation of 3.0 millimeters would break 40 percent of the fibers. It can 
be expected that this is the effective strain level desired to verify the preload 
mechanism. Because of the fraying problem associated with this method, the 
preload level was reduced to a level below which fraying occurred. After a 
series of samples were tested at various strains, it was determined that 2.25 
millimeters provided the best results for a range of room temperature curing 
times of 19 to 22 hours. At 2.25 millimeters, it can be predicted that 13 
percent of the fibers were deliberately broken. 

1. Procedure 

The finalized preload procedure consisted of : 

- Homogeneous mixing of epoxy and curing agent 

- Impregnate graphite tow with epoxy in a flat pan 

- Partially cure epoxy matrix at room temperature for 19 hours 

- Preload samples, under displacement control to 2.25mm ( .9% 
strain) 

- Final cure samples at 60°C for 16 hours 

- Tensile test samples at crosshead displacement of 10 mm/min. 

Ultimate tensile test results for the preloaded and non-preloaded 
samples are listed in Appendix D. Table XV. 

2. Interpretation 

The tensile strengths associated with the preloaded sample are 
compared with the bench mark non-preloaded samples. The non-preloaded 
samples were prepared essentially by the same procedures listed above with 
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the exception of air curing and preloading, using an automated filament 
winding set-up with no manual handling. The expected quality of the non- 
preloaded samples is comparable to well controlled production quality. The 
preloaded samples were prepared individually with all manual handling. The 
quality was expected to be lower than a fully developed production process. 
One can expect that any improvement observed with the preloaded sample 
may be realized in a fully developed process. 

The tensile strengths of both the non-preloaded and the preloaded 
samples exhibited statistical scatter precluding a deterministic comparison. 
Comparisons are made first at the non-parametric level and then at the 
parametric level together with consequences in structural efficiency and 
reliability. 

a. Non-Parametric Comparison 

Distribution-free properties of the two sets of data are 
presented for an objective comparisons. However, because of sample size 
(<25) is not sufficiently large, conclusions should be considered as 
qualitatively valid. 

The histograms of tensile strength for the non-preloaded and 
the preload sample are presented in Figure 13. Preloading process is seen to 
have no influence on the central tendency of the strength data but has a 
preceptable affect of removing the lower weak tail. The reduction of the 
lower tail is beneficial to both the structural efficiency and the structural 
reliability. Several relevant distribution-free properties are listed in Table 
IV. 
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TABLE IV. COMPOSITE STRENGTH PROPERTIES 





Mean 


Std Dev 


Medium 


Mode 


Skewness 


No Preload 


46.942 


3.07 


47.220 


47.776 


-.27 


Preload 


46.418 


2.12 


46.420 


46.424 


-.003 


Change 


-1% 


-31% 


-2% 


-2% 


-100% 



Consistent with the visual comparison, column one suggests 
no significant change of the mean strength between the two processes. 
However column five indicates that the nonpreload strength skewed to the 
weak tail whereas the skew is drastically reduced after preloading. This was 
also consistent with visual identification of lower tail shift after preloading, 
see Figures 13 and 14. 

b. Parametric Comparison 

The local load sharing model, [Ref. 3] provides the theoretical 
bases for the Weibull model for the non-preloaded samples. 
Characterization of the non-preloaded strength data by a two parameter 
Weibull model is therefore proper. Similar justification for the preload 
sample is not yet firmly established; this data can also be fitted by a two 
parameter Weibull model for expediency for comparison, with the 
understanding that the procedure is ad hoc. The Weibull parameters 
obtained through maximum likelihood estimator for the two sets of data are 
tabulated in Table V. 
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Figure 13. Histogram of Composite Failure Load Results 
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Figure 14. Weibull Distribution Plot of Composite Failure Load 



45 



LOAD (KG) 



TABLE V. COMPOSITE WEIBULL PARAMETERS 





shape 


Scale 


Nopreload 


17.48 


48.34 


Preload 


23.45 


47.43 


Change 


+34% 


-2% 



The Weibull representations are also consistent with the non- 
parametric observations, namely no significant central tendency shift as 
reflected to almost identical p values; but a significant decrease in scatter as 
reflected in the increase in a. The increase in a affects structural efficiency ( 
at a given reliability level). The magnitude of the potential improvement is 
illustrated in Table VI. 

TABLE VI. SUMMARY OF RELIABILITY AS A 
FUNCTION OF APPLIED LOAD (KGS) 



Reliability 


Nonpreload 


Preloaded 


Change 


1x10-3 


32.56 


35.33 


8.5 % 


lxlO' 5 


25.02 


29.03 


16% 


1x10-6 


21.93 


26.31 


20% 



For composite applications where functional reliability of 
1x10*3 is required, structural efficiency goes up by 8.5 %. For applications 
where national security reliability is lxlO' 5 , structural efficiency goes up by 
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16%. For composite applications where man-safe reliability ( lxl0~6 ) is 
required, structural reliability goes up by a significant 20%. 

Alternately, an increase in a can be interpreted as to cause an 
increased reliability at a given design stress level, (defined as a fraction of the 
mean strength). For example, The probability of failure (F) at 80% of the 
mean decreases from 1.2% for the nonpreloaded composite to .32% for a 
preloaded composite, a delta of 75% 

Finally, it was pointed out, earlier, that modeling the 
preloaded data by the two parameter Weibull model was ad hoc. In fact, the 
preloaded data when presented alone suggest a bimodel trend (Figure 15). 
Fitting a shape parameter , a to the lower tail return a value of a = 270 . This 
large value of a effectively provides for a structure that would be practically 
100 % reliable at any load level below the transition between the modes 
(approximately 43 kilograms) or, a 100% increase in strength at lxlO -6 
reliability. Thus potentially, the improvement from preloading is even 
greater. However, it must be understood that the number of samples tested 
herein is insufficient to establish a bimodel system. This latter interpretation 
is subject to confirmation by more extensive exploration. 
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V. CONCLUSIONS AND RECOMMENDATIONS 



As a result of individual fiber testing it can be concluded that Individual 
fiber Test system developed at the Naval Postgraduate School, NPSIT, is a 
viable testing method which produces accurate and reproducible results. 
The NPSIT test, is a one step test system which provides for the 
determination of fiber diameter and ultimate failure load. This test method 
minimizes the handling of the sample, thus providing a truer characterization 
of the fiber. An additional advantage is that the test was performed in 40 
percent less time as compared to using the INSTRON tester in conjunction 
with the Scanning Electron Microscope. 

This study clearly demonstrated that preloading a composite, prior to 
complete polymerization, has a profound effect on the lower tail of the 
Weibull distribution, hence the overall reliability of the composite. In 
addition, the optimum preload level can be independently mathematically 
estimated, based on fiber tests results. 

Further studies are recommended in the areas of; 

Use of the Micron Eye in determining fiber diameter. 

Alternate composite gripping methods. 

The effects of preloading as a function of strain rate, air 
curing time and different strains. 

Preloading composite strands at elevated temperatures. 

Statistical modeling to account for a multi -model distribution 
resulting from preload. 
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APPENDIX A. SINGLE FIBER TEST PROCEDURES 



1. FIBER SAMPLE PREPARATION 

The individual fibers were mounted on cardboard forms for ease of 
storage and handling. The samples were prepared in accordance with ASTM 
procedure D-3379-75 (reapproved 1982) [Ref 11.]. 



Fabricated cardboard mounts. (3" X 5/8"). 

Remove 10 inch strand of graphite form roll. 

Using a magnifier, carefully separate strand 
to expose individual fibers. 

Pick up one end of fiber with adhesive tape and cut to 
approximately twice the desired length. 

Lay fiber over cardboard mount and center ( Figure 16). 
Affix the ends of the fiber with tape. 

Place cardboard on storage rack and label. 

Apply Duco cement at the ends of the fiber. 



Adhesive 




Gage Length = 5 cm — ► 



Figure 16. Fiber Mount 
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2. CALIBRATION OF FIBER DIAMETER SYSTEM 



Before the fiber diameter can be determined using the diffracted laser 
beam associated with the Naval Postgraduate School Integrated Test, NPSIT 
a calibration curve must be established . 

The calibration curve used to determine the fiber diameter shows the 
results when comparing the fiber diameter as measured using the diffracted 
laser beam (NPSIT) compared to the diameter recorded when using the 
Scanning Electron Microscope (SEM). Six test specimens were mounted on a 
15 centimeter cardboard mount and then samples from each end were 
removed and carefully marked. These samples were then put under the 
Scanning Electron Microscope (SEM) and several diameter readings were 
recorded and averaged. The Scanning Electron Microscope results were 
used as a standard for determining the fiber diameter. 

Figure 17 is a graphical representation of the calibration result, with the 

equation of the straight line being. 

Diam = .89 * (NPSIT) + .3614 ( 1) 

A linear curve fit proves most useful and the results for an exponential 
or polynomial curve fit are essentially the same. 

It should be noted that as the fiber diameter increases, the difference 
between the NPSIT results and the corresponding SEM results also increases. 
By reviewing Figure 6, this can be predicted because as the fiber diameter 
increases the spacing between the nodes (X) decreases, and so does the width 
of the associated dark areas. Therefore it becomes much more difficult (and 
less accurate) to find the center of the low intensity node using the 
photocells. 



51 




3. FIBER DIAMETER MEASUREMENT PROCEDURES 

The fiber sample is first mounted in the test by inserting between the 
upper and lower clamps Then the entire stand is adjusted longitudinally to 
achieve the desire distance (L) from the photocells, (note: the spread of the 
light pattern increases as L increase). The photocell pedestals are then 
adjusted to an approximate position corresponding to the predicted distance 
between the nodes, ( ensure there is sufficient travel in the micrometers to 
allow for left and right movement) . After the pedestals are secured in place 
measure the distance between the photocells and record as C. Mount the 
micrometer in place on the photocell pedestals and zero the reading. Now the 
final distance between the photocells can be easily determined by adding C to 
the micrometer reading Delta C. 
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The final steps for determining the fiber diameter are performed with 
the lab lights off, therefore complete familiarization of all test equipment 
adjustments is mandatory. After the laser is turned on and the shutter 
opened, a red light (beam) will appear (avoid looking directly into the laser 
beam). It may be necessary to manually adjust the laser position and the 
transverse position of the test stand to obtain the optimum diffraction pattern 
at the point of interest on the fiber. Then manually adjust the height of the 
photocells to align the light pattern to the center of the photocells. The 
preliminary adjustments are now complete. 

To determine the distance between the second nodes, each photocell 
pedestal is adjusted using the micrometers. The final adjustments are made 
transversely to locate the position associated with the highest resistance 
reading on the multimeter. Readjust both pedestal separately several times to 
ensure maximum accuracy. After you secure the laser power, record the 
digital micrometer reading as Delta C. 

Calculate the fiber diameter, a, from equations (2) and (3) of [Ref 9.]. 



The accuracy of the fiber diameter is dependent upon the operators 
techniques, so five or six readings should be recorded and averaged. 



X = C + Delta C, 
a = n (lambda) L / (X/2), 
n = corresponding node (eg n = 2), 
L = 598.6 mm. C = 115 mm. 



( 2 ) 

(3) 
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4. LOAD TESTING 



Fiber failure load was determined using either the Instron Universal 
Testing Instrument Model 4206 (INSTRON) or the NPSIT system. A 
Sensotec Inc. 50 gram Load Cell Model 3 1/1435-01 was used to sense fiber 
load. The Machine Complaince test, a cross head speed of .1 millimeters per 
minute was best and for the Fiber Test a speed of 5 millimeters per minute 
(.lx gage length) was used. 

INSTRON testing is a semi-automatic process when used in conjunction 
with the HP - 85 computer and an HP 3497 A data acquisition system and the 
Interactive Data Acquisition software outlined in Appendix B. 

NPSIT testing also used the same data acquisition system , however 
manual operation of the drive motor is required. The Control Technics 
Corp. Linear Actuator Model CTC-160 (gear head ratio 262:1) and 
associated Micrometer Controller Model MMC-10 provide the same 
functions as the Instron cross head. The Speed Enter was set at -200 on the 
thumb wheel to achieve 5 mm/min and the Ratio Enter was set at 40 for a 
total travel of 45 turns in one minute. 

After load cell calibration and system compliance tests are complete, 
Fiber testing was conducted as follows. 

Run FTST program. 

Carefully mount sample in clamps. 

- Select DIAM (NPSIT only). 

Measure Fiber Diameter (NPSIT only). 

Select INPUTS and input initialize data. 
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Select ACQUIR, then carefully cut the cardboard mount 
near the bottom. 

Remove slack in fiber (NPSIT only) 

Hit program CONT and controller R/S simultaneously 
After fiber breaks and data acquisition stops, stop 
motor controller. 

Graph and store data. 

With the ENSTRON test system fiber diameter was determine from the 
NPSIT first. Also this system has automatic fiber slack removal functions. 
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APPENDIX B. INTERACTIVE DATA ACQUISITION SOFTWARE 



The Interactive Data Acquisition Software (IDAS) was designed for use 
with the HP-85 microcomputer and the HP-3497A Data Acquisition System. 
This software consists of a set of programs used in determining the ultimate 
tensile strength and diameter of an individual fiber or filament. With the 
IDAS it is possible to create a load cell calibration curve and to obtain output 
in units of grams-force. A program was also created to determine the system 
compliance in accordance with ASTM procedures. Provisions for complete 
data analysis are provided through data storage, retrieval and plotting 
routines. The software gathers data at a rate of 600 points per minute and 
records and plots Load vs Displacement of the fiber sample. 

The IDAS can be used with the INSTRON Model 4206 materials tester 
or the Naval Postgraduate School Integrated Test system (NPSIT). The 
programs listed are associated with the INSTRON test system. When using 
the NPSIT system, the program steps that provide for automatic control of 
the INSTRON (OUTPUT "704") must be deleted. Crosshead control and 
pretensioning of the fiber is accomplished manually. A program sub routine 
is listed in the FTST module and is used to determine fiber diameter using the 
NPSIT system. Functions keys are defined to provide ready access to 
commonly used sub routines. 

Mr. Jim Nageotte served as the principle computer technician during the 
development of these programs. 

Complete listing of the IDAS programs and a brief description follows. 
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C 3 = S 1 <K'.J/ r X5 



OCT 26 



" h i ri •• 
h < i :• 

"Hi 1 M " 
B C T ) 



’’LDCALB” calibrates the load cell to 
convert the voltage out into grams 



ON " 

ON 

Select the number of calibration points. ( six 
is optimum, maximum ten) 

Start of loop, once through for each data 
point. 

' 3 V 0 Input load level applied to load cell. 



The 3497 A reads the output voltage A(I). 
And excitation voltage B(I). 



Fitting the data to a linear function. 
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520 C3= 100*C3 



530 


C3=RBS ( C3> 




540 


PRINT 


" L 0 R D = 


" ; L < K ) ; " ?m" 


550 


PRINT 


"MERN 


V = " > X 5 < K ) ; " 


560 


PRINT 


" S 0 


SI <k> " V 


570 


PRINT 


" f: i 


• = " ; C3 


5 SO 


PRINT 


" L C . f 


IJTPIJT VOLTflGF 


5 9 0 


PRINT 


"MERN 




6 0 0 


PRINT 


"S D 


»J — •• ; C; Ci , '« l (( l " 


610 


PRINT 


" C ' 


V > = *' ; r 9 


620 


PR I NT 


"EXCITATION uoltrge 


630 


PRINT 


m MEhN 


| | _ M . £ O . II |J II 


640 


PRINT 


" s . D . 


V=" ;S3;" V" 


650 


PRINT 


" C U 1 


■;> = " ; cs 


66 u 


PRINT 


It II 




6 70 


NEXT K 






6 S O 


X 1 = 0 0 


Y 1=0 




690 


X 3 = 0 i? 


X 2 = 0 




7 O O 


FOP K = 


1 TO Ml 


7 10 


X1=X 1 +UOO 




720 


Y 1 =Y 1 +• L < K> 




730 


X 3 = X 3 + V < K 


< k ;■ 


740 


X2=X2+ 


vao-2 





Begins printing out statistical parameters. 



Continues once for each data point (load 
level). 



750 


NEXT K 




76 y 


X 1 =X 1 '-'N 1 




7 70 


X 2 = X 2 N 1 




7S0 


X 3 = X 3 s N 1 




790 


Y 1 = Y 1 •• N 1 




800 


D 1 = X 2 - X 1 2 




S 1 O 


r 1 = < X 3 - X 1 * Y 1 :t 


D 1 


820 


C 2 = ■; X 2 * Y 1 - X 3 * X 


1 >^D1 


830 


PRINT "L=R*V+B 


li 


840 


PRINT " R = " f. 1 




850 


PRINT " B = " C 2 




860 


R6 = 0 




8 70 


FOR K = 1 TO N 1 




880 


G2=<L(K>-C1S V< 


K > - C 2 



For least squares method, voltage = V(K) and 
Load = L(K) 



After all load levels are recorded, the program 
calculates and returns the best fit 
coefficients, A and B. 



Determines the goodness of fit which is 
expressed by the "residual". 



8 9 O R 6 = P 6 + G 2 
900 NEXT K 

910 PRINT " RES I DURL = " , RS The graph routine plots the load in grams 

920 DISP •'SCREEN GRfiPH ROUTINE" versus voltage out of the load cell. This 
930 DISP " M R X . L 0 R D = " graph appears on the computer screen and is 

940 I N P U T Y 2 printed by the thermal printer. 

950 DISP "MIN. LORO=" 

960 INPUT VI 
970 U2=X5< 1 } 

980 U1=X5<1> 

990 FOP 1=1 TO N1 

1000 IF X5 < I > > U2 THEN U2 = X5 (. I ) 

1O10 IF X5 < I > <U 1 THEN U1=X5(I> 

1O20 NEXT I 

1O30 PRINT "X MRX U 000 V/V > = " j U2 

* 1 GOO 

1040 PR I NT " X tl IN <1 000 V/U > = " ; IJ 1 

* 1 0 O O 

1050 DISP "INPUT MRX FOR X=" 
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10*0 INPUT U2 

1070 CUSP "INPUT MIN FOR X=" 

1G80 INPUT U1 

1090 FOR 1=1 TO Ml 

1 1 00 B 1 * I > = X5 < I > *. 1 000 

1110 B2< I :>=S1 < I > t 1000 

1120 NEXT I 

1130 81 = 160 < U 2 - U 1 > 

1140 S2= ! 40 ' < Y2-Y l > 

1150 FOP 1=1 TO 5 
11* 0 tj F B P H 0 ij L LEAP 
1170 LDIP O 

1 1 80 SCALE -48 . 208 . -7*, l SO 
1190 X fl X I S 0 .• 3 2 O , 1 * 0 
1 200 VAX IS 0.28 0.1 40 
1210 PEN 1 
1220 PE HUP 
1 2 30 W 1 = 1 * 0 ' < IJ2-U1 > 

1240 W 2 = 1 4 O / ( V 2 - V 1 
1250 FOR 1=0 TO 5 
12*0 U = U 1 + ( IJ2-U 1 ,• r 5 1 1 
1270 P3='U-U1 > fcHl 
1280 MOVE P3,-10 
1290 LABEL V ALT 1 U,' 

1300 NEXT I 
1310 FOR 1=0 TO 5 
1320 Y = V 1 + < Y2-Y 1 V SU 
1 330 03= < Y-Y 1 ^ tW2 
1340 M 0 W E -23 -Q? 

1350 LABEL UAL# ( Y ) 

13*0 NEXT I 

13 70 MOVE 4 0.. -2-: 

1380 LABEL " lOOOtVout/Vexc " 

1390 DEG 

1400 MOVE -29.25 

14 10 L D I F' 90 

1420 LABEL "FORCE (•?!») " 

1430 LDIP 0 

1440 FOR 1=1 TO N 1 

1450 Z 1 < I .> = < B 1 < I > -U 1 ' *W 1 

1460 Z2< I ) = < L I >-V 1 >*W2-2 

14 70 MOVE 2HI)..Z2(I) 

1480 LABEL 

1490 NEXT I 

1500 G 1 = C 1 / 1 O 0 0 

1510 FOR I=0 TO 101 

1 520 U=U 1 + •' IJ2-IJ 1 > •" 1 DO* I 

1530 Y=G1*U+C2 

1540 IF Y < Y 1 OR Y>Y2 THEN 1580 
1550 Zl(I) = <U-U!mil 
15*0 Z2< I ') = <Y-Y1 >*M2 
1570 PLOT Z1 •' I >, Z2< I > 

1580 NEXT I This rouline plots the graph on the HP 

1590 C 0 P Y Plotter. 

1*00 DISP "00 YOU WANT A HARD CO 
PY ON PLOTTERS <Y/N>" 
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1610 
1620 
1 630 
1640 
1650 

1 660 



INPUT KS* 

1F K8t=«N" THEN 2170 
PRINTER IS 10 
CONTROL 10-5 ; 4ft 
OUTPUT 10 "I H ; SP 1 ; I P2400 , 1 
600 .• 8S00 , 6800 : " 

OUTPUT 10 • "SCO- 1000.0 • 1000 



1670 OUTPUT 10 " PUO • 0PDO • 1 000 > 1 

000 , 1 00O .• 1 000 ■ 0 • 0 •- 0PU " 

1 680 W 3 = 1 0 O 0 '' U 2-U 1 > 

1680 W4=lO00 y < Y2-Y1 > 

170O OUTPUT 10 "810 . 2.- 0 3; TL 1 5 

. 0" 

1710 FOR 1=0 TO 5 
1720 Y=Yl+< Y2-Y1 > '541 
1730 Y4 = < Y-Y 1 > 4W4 
1740 y=INTCY> 

1750 v 4 = I N T ( Y 4 > 

1760 OUTPUT 10 ; "PR 0 , " , Y4 , " YT , " 
1770 OUTPUT 1 0 ; " CP-5 , -8 . 07 ; LB " 

Y ; CHP$(3 ) 

1780 NEXT I 
1700 FOR 1=0 TO 5 
1 880 U =IJ 1 + < IJ 2-U 1 > 4 I s5 
1810 U4=vU-Ul >4W3 
1820 U4 = IHTMI4> 

1830 U=INT<U> 

1840 OUTPUT 10 ; "Pfl"i U4, “ ,0; XT; " 
1850 OUTPUT 10 "CP-1 .3,-1 ; LB" IJ 

CHR*<3> 

I860 NEXT I 

1870 OUTPUT 10 ; "SI .30, .42" 

1880 OUTPUT 10 ; "FR400 , 0 ; CP-2 > -2 
.3; LB 1 0004Vou t /Vexc" ; CHR*<3 

> 

1 830 OUTPUT 1 0 " PRO , 460 ; D 1 0 , 1 ; C 

P- 2. 6, 2. 6; LB FORCE (<?mp ";C- 
H R f ( 3 > 

1 800 OUTPUT 1 0 " 0 I ; PU " 

1810 FOR 1=1 TO Ml 
1820 Z 1 < I > = ( B 1 <. I > -IJ 1 > 4W3 
1330 Z2< I :> = <L< I >-Y 1 >4W4 
1840 ZH I >=IMT<21 <1 > > 

1850 Z2< I >=INT<Z2< I >> 

1 360 OUTPUT 1 0 ; " F'U , " ; " Pfi " , Z 1 < I > 
, Z2( I ) ; "CP-0 .1,-0 . 2; LB*" ; CH 
R$<3> 

1370 NEXT I 

1980 OUTPUT 10 ;"PU" 

1990 FOR I=0 TO 101 
2080 1 1= IJ 1 + < IJ2-IJ 1 10041 

2010 Y=G 1 4U+C2 

2020 IF Y < Y 1 OR Y > Y 2 THEN 2060 
2038 Zl( I >=<U-U1 >4U3 
2040 Z 2 < I > = < Y - Y 1 > 4 N 4 
2850 OUTPUT 10 ; "Pfl" ; Z 1 < I > , Z2 < I ) 
; "PD" 
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ro ro ro ro ro ro r ro ro rv ro ro 



060 NEXT I 

070 OUTPUT 10 "PU0, 900, 100- 900 

ii 



0 S O 
090 

1 00 

1 10 
1 20 
130 



OUTPUT 10 M P U " 

OUTPUT 10 ,"31.22, 38" 

DISP "ENTER THE LEGEND ENT 
ER ' 0 ' TO EXIT" 

INPUT P?$ 

IF F' 7 $ = " 0 " THEN 2150 
OUTPUT 1 0 : " CP LB " .; P7$ ; CHRJ 



The user can enter a label or legend on the 
finished plot. 



140 L> U T 0 2 1 O 0 
150 PRINTER IS 
160 DISP " 

170 END 



"LDCALB" ends. 

END LDCflLB " 
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2. MACHINE SYSTEM COMPLIANCE MODULE (CTST) 



10 ! CTST 8 7 

20 ! JIM NAGEOTTE -JAN 13 37 

3 0 1 

40 CLF-fiR 

50 0 1 3 P <? 0 1 3 P •? CUSP !? DISP (5 
CUSP " CTST87" 

60 CUSP e CUSP " USE THE FUHCT 
ION KEYS." 

70 PP INTER IS 2 

SO 0 I M P 1 450 > , L <450 1 , D <450 > , Z 1 < 
4 50 > .22 1 4 50 > 

SO PRINT “IF PROGRAM HALTS, TYF 
E > ' CONT 100’ TO RETURN TO ME 
NU . " 

100 PRINT 
110 GOTO 130 
120 CLEAR 



1 30 


ON 


KEY# 


4, 


“END" GOTO 220 


1 40 


ON 


KEY# 


1 , 


" INPUTS” 


GOSUE 240 


150 


ON 


KEY# 




” ACQU I R" 


GO SUB 530 


160 


ON 


KEY# 


• 1 ' f 


"CUSP" GO 


SUB 1110 


1 70 


ON 


KEY# 


6 ■ 


'•RD_DRT H 


GOSUE 156 


ISO 


ON 


KEY# 


i 


” FLOTTP” 


GOSUE 186 


1 80 


yj 

ON 


KEY# 


5 , 


"SAU.DAT” 


GOSUB 17 



1 0 

2O0 KEY LABEL 
210 GOTO 130 

220 CLEAR 0 CUSP "END CTSTS7" <3 
END 
2 30 ! 

240 ! 44 INPUTS 44. 

250 CLEAR 

260 1 0 ISP” ENTER DATA RATE < p t 
l n > ” 

270 ! INPUT P 

280 R=6O0 ! PRESET DATA RATE 

280 DISP “ENTER ELAPSED TIME <30 
OP 4 5 > " 

30O INPUT E0 

310 CLEAR 8 DISF "ENTER SPEED IN 
MM. "MIN ” 

320 INPUT S 

325 PRINT "SPEED = " ; S ; " NM/"M I N " 
330 E 1 = E 0 4 1 0 0 0 

340 DISP "ENTER DATA FILE NAME" 
350 INPUT NT 
360 CLEAR 

376 CUSP "P-A4V+B" S DISP 
388 DISP "ENTER ’A’ FROM LOCALE 
LISTING" 

382 INPUT C3 

400 DISP S DISP "ENTER ’E’ FROM 
LOCALE " 

410 INPUT E 



The system Compliance module or "CTST" 
measures the compliance or displacement of 
the system without a fiber sample. 



This program is menu driven using function 
keys to access the sub routines. Should 
execution halt, return to the menu by 
typing "cont 130". 



The keys are defined here. 



The program will continue to loop on the 
menu until a function key is depressed. 



The "INPUTS" sub routine is required first. 
Input the variables and information needed 
along with the values of A and B obtained 
from "LDCALB". 
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4 20 

439 

440 

450 

4 1? 0 

4 7u 

4 8 0 

4 90 

500 

505 

510 

520 

5 3 O 

540 

550 

560 

570 

590 

600 

6 1 0 

620 

636 

640 

650 

653 

655 

658 

6 6 0 

661 

662 

6 6 3 

664 

665 

666 

667 

670 

6 75 

680 

685 

687 

688 

69 1 

6 9 4 

695 

700 

710 

720 

730 

740 

750 

760 

765 

770 



CLEftR 709 

OUTPUT 709 "ft! 10" 

ENTER 709 .; V 
OISP " EXC I T ft T I ON UQLThGE 

: V 

ft=ftBS •: 3 * 

ft= I NT ( hT. 1 000 > 1 000 

01SP "ft=" ft 
DISP " B = " • B 
COPV 
CLEftP 
RETURN 

i 

■ 4 4 ft COM I F 44 

IF E0>© THEN 590 
BEEP 8 BEEP ft BEEP 0 CLEftR 
DISP "’INPUTS' MUST BE RUN F 
I PST i " 

RETURN 

CLEftP 

BEEP 

FOR 1=1 TO 7 ft DISP ft NEXT : 
DISP "CHECK IEEE BUTTON- MOL 
NT SftMPLE FIND BURN " 0 OISP 
DISP "HIT CONT WHEN RERDY" 
PftUSE 

CLEftF: O DISF 

DISP "FINDING ZERO LOF>D" 
OUTPUT 704 .: " K 1 3 ; S 

OUTPUT 709 i "ftRVRl VT3VO5VFl0> 
E2 " 

OUTPUT 709 ."ft 10" 

ENTER 709 .; L 
F' = ft 4 L + B 

IF P<9 THEN GOTO 670 
OUTPUT 704 ; " K 2 " 

WRIT 500 

OUTPUT 704 "K0" 

GOTO 660 

CLEftR 0 BEEF 0 DISP 0 DISP 1 
MOVING UF " 

OUTPUT 704 . "K3" 

OUTPUT 709 " ft 1 0 " 

ENTER 799 .: L 
P = fl T L + B 

IF P < . 06 THEN GOTO 680 
OUTPUT 704 :"K21" 

BEEP 0 DISP "TIMER ON" 

ON TIMER# 1 / E 1 GOTO 760 
1 = 1 

OUTPUT 709 
ENTER 709 .; L<I) 

1 = 1 + 1 
WRIT 23 
GOTO 710 
OFF TIMER# 1 
OUTPUT 704 . " K 0 " 

1 = 1-1 



The HP 3497 A records the excitation 
voltage applied to the load cell. 

S " The equation that translates the load cell 
voltage into physical units of grams is 
adjusted for the present excitation voltage. 

"INPUTS" is terminated by returning to the 
menu. 



The data acquisition sub routine is title 
"ACQUIR". Acquir checks die time duration 
selected and if it finds it to be not greater 
than zero, prompts the operator to run 
"INPUTS" and returns to the menu. 



' The program pauses and waits for the 
1 operator to hit "cont". At that time a timer 
is started and data is acquired at a rate of 600 
data points per minute. 



Note: OUTPUT ”704" provides for 

automatic control functions to the 
LNSTRON tester. These functions provide 
for a slight tension in the fiber prior to the 
start of data acquisition. 

The start of the timer is defined at the zero 
load. 
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7 S 0 WRIT 50 
790 CLERR 

£:00 FOR T= 1 TO 7 i? 01 SR G NEXT T 
S10 BEEF' 

830 01 BP 

3 4 0 clear G DISF " WRIT" 

350 H= I 

860 PRINT G PRINT H j "DATA POINTS 

II 



8b 5 PRINT 

3b b OUTPUT 704 ; " K I " 

870 FOR 1=1 TO N 
S30 P < I ) = Af L < I > + B 
398 NEXT I 
910 P2 = Pc:i> 

920 P 1 =P< 1 > 

930 FOR 1=2 TO N 
940 IF F < I > > P 2 THEN P2 = P<n 
950 IF PMXF'l T HEN P1=P<I> 
970 NEXT I 

980 PRINT "p- MAX :■ = " .; P 2 
990 PRINT " P C M 1 N > = " ; P 1 
1090 PRINT 
1030 CLEAR 
109O RETURN 
110O ! 

1110 ! tt OISP GRAPH ON CRT 
1120 PRINT " ";Nt 

PRINT 



The voltage is converted to grams force 

(P(I)). 



"ACQUIR” terminates by returning to the 
menu. 

44' 

G 



1130 CLERR G OISP "INPUT P<MAX> 

FOR V A X I 8 1 The function key defined as "DISP” calls the 

1 140 INPUT P2 sub routine which plots the data on the 

1150 P 1 =0 screen and on the thermal printer. 

1160 GCLEAR 
1170 GRAPH 
1180 LDIP O 

1190 SCALE — 4 8 .• 2 0 8 ■ — 3 b , 1 5 b 
1 20O X A X IS 8,32.0-160 
1210 VAX IS 0,23,8,140 
1220 PEN 1 
1238 PENUP 

12 48 31 = 1 bO.' <. E O / 60 > 'S 

1258 S 2 = 1 4 0 P 2 

1268 FOR T= 1 TO 5 

1278 X = S 4: T / 1 O 

1280 K= I NT C X4 1 080 ) , 1 000 

1290 X 2 = X 4 S 1 



1300 MOVE X2.--10 
1310 LABEL VAL*<X> 
1328 NEXT T 
1338 FOR T=0 TO 5 
1340 P=P2/5*T 
1350 P3=P4S2 
1360 MOVE -23- P3 
1370 LABEL VALSHPT 
1330 NEXT T 
1390 MOVE 40,-23 
1400 LABEL "DISP Onrn> 
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14H3 
1420 
1430 
1440 
1450 
1 460 
1470 
1480 
1 400 
1500 
1510 
1520 
1 5 30 
1540 
1 550 
1 580 
1570 
1 580 
1 500 
1 6 0 0 
161 0 
1 620 
1630 
1 640 
1650 
1 6 6 0 
1670 
1630 
1600 
1 700 
1710 
1 720 
1 730 
1 740 
1 750 
1760 
1 770 
1780 
1 700 
1 800 
1 8 1 0 
1820 
1 8 3 0 
1840 
1850 
I860 
1870 
1880 
1 800 
1800 
1010 
1 020 
1030 
1040 
1950 



DEG 

HOVE -29.25 
LDIP 90 
LABEL " FORCE 
LDIP 0 

FOR 1=1 TO H 
D ( I >=S*< I /P ' 1 
zi(n=D(i) •! s l 
Z2C I >=P< I > *82 
PLOT Z1<I).Z2U> 

NEXT I 
C 0 P V 
CLEAR 
RETURN 

i 

! ft READ DATA ** 
CLEAR 

CUSP "INSERT DATA TAPE" 
CUSP "HIT CUNT" 

PAUSE 

ASSIGN# 1 TO Hf 
READ# 1 : t!>EO»S 

FOR 1=1 TO N 
READ# 1 s PU) 

NEXT I 

ASSIGN# 1 TO :t: 

CUSP “ DATA READ" 

R = 6 O O 

BEEP 

RETURN 

! tt 8 A 1 1 _ D A T ** 

CLEAR 

IF P < 3 > <" > O THEN GOTO 1700 
BEEP 15 BEEP 6 CLEAR 
CUSP "NO DATA IN MEMORY" 
CUSP "HIT CONI" 

PAUSE 

RETURN 

DISP "INSERT DATA TAFE" 
CUSP 

DISP "HIT CONT" 

PAUSE 

CLEAR 8 DISP "WAIT" 

M = ( H + 1 >4 8 + 50 
H = I N T ( M / 2 5 6 > + 1 
CREATE H*,H,256 
ASSIGN# 1 tq 
PRINT# 1 ; N,E0,S 

FOR 1=1 TO N 
PRINT# 1 ; P < I > 

NEXT I 

ASSIGN# 1 TO 4 
CUSP "DATA SAVED" 

RETURN 

i 



D(I) is ihe displacement 



"DISP” terminates by returning to the 
menu. 



The load data is read in. 

The number of data points saved, the time 
the test ran and the speed are inputted.. 



A prompt and beep indicate the data has been 
read in. 

The sub routine is terminated. 

The "SAV-DAT” routine is just the counter 
part of the above routine. It first checks to 
see that there is data in the memory to be 
saved. 



The memory space is computed and the file 
opened. The data is then saved. It is in the 
same form as that read in by "RD-DAT". 



The sub routine is terminated 
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** PLOTTER *:* 



I 960 ! 

1970 CLEhR 
1930 DISP "ENTER V MRX FOR Y 
s " 



The "PLOTTER" function is design to create 
P! j a plot on the HP Plotter. 



1990 INPUT Y 2 



2 0 0 0 


DISP 


"ENTER Y MIN 


FOR 


v HXI _ 


20 1 0 


C N 

INPUT 


Y 1 


Select the maximum and minimum values 
for the X and Y axis on the plot. 


202 O 


DISP 

II 


"ENTER X MRX 


FOR 


X HXI 


2030 

2040 


INPUT 

DISP 


V ? 

"ENTER X MIN 


FOR 


X HXI 


2050 


n 

INPUT 


XI 




This plot will be on the HP plotter similar 
to the plot on the CRT screen. 



2 0 6 0 S 3 = 10 0 0 / ( X 2 - X 1 > 

2 0 70 S 4 = 1 U 0 O / (. V 2 - Y 1 > 

2O30 PRINTER IS 10 
2090 CONTROL 10-5 ; 43 

2 1 00 OUTPIJ T IQ ; " I N : SP 1 ; I P24 OO , 1 
6 0 0 .. 0 3 0 0 . 6900 : " 

2110 OUTPUT 10 " SCO , 1 ©CO ,0,1 003 

it 



2 1 20 

2130 

2140 

2150 

2130 

2170 

2130 

2 1 9 O 
2200 
2210 
2220 
2230 
2240 

2250 
2 2 6 0 
2270 



2 2 3 0 



2290 

23O0 

2310 

2320 

2330 

2340 

2350 

2360 



OUTPUT 10 ; " F'lJO , OP DO • 1 00O , 1 

000 , 1 000 .• 1 000 • 0 , 0 , CPU " 

OUTPUT 10 ; "SI0.2,0.3;TL1 5 
, O " 

FOR 1=0 TO 5 
Y = I N T < Y 1 + < Y 2 - 7 1 >/5*I ) 

Y4=INT< (Y-vi yi: S4> 

OUTPUT 10 ; " PR 0 , " , Y4 , " YT .; • 

OUTPUT 10 "CP-5, -0 . 07.: LB" j 

Y ; CHRf < 3 > 

NEXT I 

FOR 1=1 TO 5 
X=X 1 + < X2-X 1 >/5*I 
X 4 = I N T O. X - X 1 >*S3> 

OUTPUT 1 O ; " PR " ; X4 , " , O :■ XT ; " 

OUTPUT 10 : "CP-1 . 3, -1 ; LB" ; X 

} CHRf -.3 > 

NEXT I 

0 IJ T P U T 10 .; "SI 0 > -4 2 11 Th e fables for the axis are put in place. 

OUTPUT 10 " PR220 , 0 ; CP4 , -2 . 

3; LBO I SF'LRCEMENT (mm) " ; CHRf 



c 3 > 

OUTPUT 10 ; " PRO • 460 D I 0 , 1 C 
P-2 . 6, 2 . 60; LB FORCE Cgrr.V'iC 
HRf < 3> 



OUTPUT 10 "DIjPU" 
FOR 1=1 TO N 
D< I ^Stl'P 



Here displacement is represented by D(l) 
while P(I) is load. 



Z 1 < I ) = I NT < < 0 < I ) -X 1 > * S 3 > 



Z2 <1 > = I NT < P < I > - Y 1 > YS4 > 
OUTPUT lO ; "PR" , Z1 < I ) , Z2< I ) 

. •• P p " 

NEXT I 

OUTPUT 1 0 , " PU0 , 900 , 1 00 , 900 
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2370 
2 .3 8 y 
2390 

2 4 O 0 

2410 

2 4 2 0 
24 30 
2440 

2450 

2460 

2470 

2535 



0 U T P IJ T 10 '* p |j '* 

OUTPUT 10 " ft 1 2" ? . 3 ft " 

DISP " LOCft T E* PEN" FOR' LEGEND 

" After the plot is completed, the pen moves 

DISP 2 D I S F "ENTER THE LEGE 101 * 10 u PP cr kit corner where a legend may 
HD YOU WAN 7 . " be created. 

DISP "ENTER ft 'O' WHEN DONE 



INPUT Pi 

IF P i = " 0 " THEN 2460 
OUTPUT 10 i " CP LB" ; Pi CHRi < 



O 0 TO 24 0 0 
PRINTER IS 2 
RETURN 
CLEAR 709 



The program is completed and returns to the 
menu. 
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3. QUADRATIC FIT MODULE (QFIT) 



1 w 



20 
30 
4 0 



50 
60 
7 0 
80 
90 
1 00 



1 10 
1 20 
1 30 
140 
150 
1 60 
170 
ISO 
1 90 



2 0 O 
210 
220 
230 

24 0 
250 
260 
270 
2 y 0 

290 

300 

310 

320 

330 



340 

350 

360 

3 70 
380 
380 
400 

4 10 
420 
4 30 
4 40 
441 



' XI Ur I T 10 -24/8 

b 

' JIM NAGEOTTE 

i 

' FIT the ML TXT DATA BY H Q 
UADRAT I C EQUATION 

i 

CLEAR 

PRINTER IS 2 

SHORT P < 450 :• , Z 1 <45© > , Z2C450 ) 
D IN E < 2 • 2 ^ > F < 2 ' G < 2 ) • D < 4 5 O ) 
BEEP 0 DISP “INSERT DATA TAP 
E" 

DISP "INPUT FILE NAME" 

INPUT NT 
ASSIGN# 1 TO NT 
READ# 1 N .• E 0 • S 

FOR 1=1 TO N 
READ# 1 ; PCI) 

NEXT I 

ASSIGN# 1 TO * 

P — 6 0 0 ! DATA SAMPLING RATE<p 
n t s / rn i n > 

FOR 1=1 TO N 
Dv I :> = r*s/p 
NEXT I 

DISP "DO YOU WANT TO DISP TH 
E GRAPH ON SCREEN (Y/N)?" 
BEEP 

INPUT AT 

IF AT = " Y " THEN GOSUB 1020 



"QFIT” fits a quadratic equation to the 
compliance data obtained from "MCTST". 



Input file name. Use the file stored by 
"MCTST". 



A graph or plot at this point would give 
results identical to the output of "MCTST". 



i 



DISP "DO YOU WANT A HARD COP 
Y OH PLOTTER? < Y/N > " 

INPUT AT 

IF A T = " Y " THEN GOSUB 1600 



DISP "Q FIT SUBROUTINE" 
DISP "COMPUTING ABOUT 2 
UTES " 

E < O , O > =H 

E < 0.1) =0 

E < 0 .• 2 ) =0 

E < 1 ■ 1 ) = 0 

E < 1 .• 2 ) =0 

E < 2 > 2 ) =0 

FOR 1=21 TO N 

E<0, 2>=E<0, 2)+F< I )-2 

NEXT I 

FOR 1=1 TO N 
E-..0, 1 >=E<0, 1 > + P< I ) 

E < 1 , 2 > =E < 1 , 2 > +P < I > -'2* < P < 

N 



MIN 

The "QFIT" sub routine requires about two 
minutes. The exact time depends on the 
number of data points to be fitted. 



Initializing the matrix 



I > t H Least square method 
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52 O E 1 2 • 2 >=E < 2 , 2 > + P ( I ) ^2* <P<I '••'•2 
i N > 

53U next i 



535 


E < 1 


, 1 


> = 


E< 


0 • 2 


> 








540 


E < 1 ■ 


1 


j rz 


Ev 


i . 2 


'l 


/N 






550 


E < 2 , 




•» - 


E < 


v , v 


l 


/N 






5 7 0 


g 1 0 ; 


> = 


fl 


A 


G < 1 


> 


- 0 1 5 


G 1 


2 > = 0 


v?0 


FOR 


I 


= 1 


TO N 










6 6 0 


G UT : 


1 = 


G< 


0 > 


+ 0 ■ 


I 


) 






630 


G < 1 :. 


« = 


G ( 


1 > 


+ DC 


1 


) :i P * 


I > 




660 


G • 2 : 


1 = 


G ( 


2 > 


+ D 1 


I 


> i P'. 


I )' 





670 NEXT I 

680 CUSP "MRTPi: OPEPflT I ON •' 

696 E< 1 • 6> = E • 6 .. 1 > 

7 0 0 E < 2 • 0 '• = E ( 0 . 2 •• 

710 E C 2 . 1> = E < 1 , 2 > 

726 MAT DISP E 
730 MHT F = SYS C E • G.' 

740 PR INTER IS 2 

750 PRINT ” ":H$ 6 PRINT 

760 PRINT ,, P = E + F*DELTR+G*0ELTR''2 

It 

770 PRINT F-. 0.) 

780 PR INT " F = " : F v 1 > 

790 PRINT " G = " ; F C 2 > 

800 CUSP "P = E + F*0ELTFi + G*0ELTfl''2" 

8 10 DISP " E = " ; F ' 0 ) 

820 DISP " F = " ; F 1 1 > 

830 DISP "G=":F<2> 

840 MRT PRINT F 
856 E = F < 0 > 



By symmetry 



The best fit coefficients E, F, and G arc 
returned and printed out. 



860 F = FM> 

870 G=F<2> 

886 CLERR @ DISP "COMPUTING" 

890 FOP 1=1 TO N 

906 D < I > =E + FtP( I > + G TPC I ) ~2 

916 NEXT 1 

920 BEEP 



The coefficients arc then used to calculate the 
estimated curve of displacement as a function 
of load P(I). 



930 CUSP "DO YOU WANT ft COPY OF 
ESTIMRTED CURVE ON THE SCREE 



940 

950 

960 



N ? •' Y / N > " 

INPUT Hf 

IF A f = " Y " THEN GOSUB 1020 
DISP "00 YOU WANT ft HRPD COP 



A plot of this curve , shows a smooth 
representation of the original curve. 



Y OF ESTIMATED CURVE ON PLOT 





TER? C Y 


••••' N > " 






9 7 0 


INPUT 


H* 






980 


IF A* = 


"Y" THEN 


GOS 


UB 16O0 


990 


BEEP 








1 000 


DISP 


"END OF 


OF IT 


a 


1010 


END 








1 020 


! tt 


SCREFN 


PLOT 


** 


1030 


DISP 


"DISPLAY 


THE 


RESULT 




THE 


SCREEN " 






104 0 


PRINT 


il 




" ; N$ 


1050 


PRINT 








1060 


Y2 = P < 1 > 







The main program ends. 



The screen plot routine is the same as that 
used in "MCTST” and plots the graph on the 
CRT screen and thermal printer. 
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1070 


V 1 = P <1 :■ 




10SC 


FOR 1=2 TO H 




1 890 


IF P < I > > Y 2 THEN Y2 


= P< I > 


1 1 00 


IF P •' I > < Y 1 THEN Y 1 


=P(T' 


1110 


NEXT I 




1 1 28 


PRINT " Y MAX=",Y2 




1 138 


DISP "INPUT Y MAX 

C* ** 


FOR Y 


1 1 48 


INPUT Y2 




1 150 


PRINT "Y H I N = " ; Y 1 




1 1 88 


DISP “INPUT Y MIN 

ii 


FOR Y 


1 170 


INPUT VI 




1 180 


DISP "ENTER X MRX 

M 


FOR X 


1 1 90 


INPUT X 2 




1 208 


DISP "ENTER X MIN 

C li 


FOR X 


1 2 1 8 


INPUT XI 




1228 


GRAFH 




1 2 3 U 


GCLERR 




1248 


LDIP 0 




1258 


SCALE -48, 288 • -36, 


156 


1 260 


XRXIS 0,32,0,160 




1270 


YRXIS 0,28. 0-14O 




1 280 


PEN 1 




1 298 


FENUP 




1 300 


X 1=0 




1318 


S 1 = 1 60 (' X 2 — X 1 1 




1320 


82=140 • ' < Y 2 - V 1 ’> 




1338 


FOR 1=1 TO 5 




1 34 8 


X=X2 y 54 I 




1 350 


X3= CX-X 1 > 4S 1 




1 368 


NODE X 3 • - 1 0 




1 370 


LABEL VflLf<X> 




1388 


NEXT I 




1 3 9 0 


FOR 1=0 TO 5 




1400 


Y=Yl+< Y2-Y1 JX5*I 




14 18 


Y3=< Y-Y 1 > 4S2 




1420 


MOVE -23, Y 3 




14 38 


LABEL VAL$ < V ') 




1440 


NEXT I 




1450 


MOVE 40.-23 




1468 


LABEL "DISPLACEMENT (nm 


1 478 


DEG 




1488 


MOVE -29,25 




1 498 


LDIP 90 




1500 


LABEL "FORC E < ?m ) " 




1510 


LDIR 0 




1520 


FOR 1=1 TO N 




1 530 


2 1 < I > = < 0 < I ) -X 1 > t-S 1 




1540 


Z2< I >=<P< I >-Yl)*S2 




1550 


PLOT Z1 < I > , Z2< I > 




1560 


NEXT I 




1570 


COPY 




1588 


RETURN 




1590 


i 




1608 


! 44 PLOTTER 


** 



FIX I 

FIX I 
FIX I 
FiX I 



The X axis or displacement is represented by 
D(I) while the Y axis or load is represented 
byP(I). 

The plotter routine plots the graph on the 
HP Plotter. 
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1610 


UISP 

V 


"ENTER Y M R X " 


£ 


INPUT 


1 620 


DISP 
Y 1 


"ENTER 


Y MIN" 




INPUT 


1 630 


OISP 


"ENTER 


X MAX" 


e 


INPUT 


1 640 


oTsp 


"ENTER 


X MIN" 


e 


INPUT 



XI 

1*50 83 = 1 0 9 0 '( X 2 - X 1 > 

1 *60 84=1 O 0 O '' C Y 2 - Y 1 > 

1670 PRINTER IS 10 
1680 CONTROL 10,5 ; 48 
1690 OUTPUT 10 " IN; SP1 ; IP2400, 1 

600 , 8800 • 69O0 ; " 

170O OUTPUT lO ; "SCO , 1 00O ,0,1 00O 

1710 OUTPUT 10 ■ " P U 0 , 0PD0 , 1 0O0 , 1 
000 , 1 000 .• 1 O0O . 0 . 0 , UPU " 

1720 OUTPUT 10 • "SI0 . 2, 0 . 3; TL1 . 5 
, O " 

1730 FOP I=0 TO 5 
1 740 Y= I NT < V 1 + < Y 2 - Y 1 > ✓ 5*1 > 

1750 V 4 = I N T < Y - V 1 ) * S 4 > 

1 760 OUTPUT 10 ; " PH O , " , Y4 . " YT ; " 
1770 OUTPUT 10 • "CP-5, -0 07; LB" ; 

Y; CHPT<3> 

1780 NEXT I 
1790 FOR 1=1 TO 5 
1 800 X = X 1 + < X 2 - X 1 } '' 5 * I 
1810 X4= < X-X 1 ) *83 
1820 X 4 = I N T >' X 4 > 

1830 OUTPUT 10 ; "PA", X4 , " , O ; XT ; " 
1840 OUTPUT 10 ; " C P - 1 3,-1; LB" ;X 
; CHR-T < 3 > 

1850 NEXT I 

I860 OUTPUT 10 "SI. 30,. 42" 

1870 OUTPUT 10 , "PA370, 0; CP-2, -2 
. 3; LBDISPLACEMENT (mm) " ; CHR 
* f. 3 "> 

1 880 OUTPUT 10 " PRO ,460-010, 1 : C 

P-2.6,2.68 , LB FORCE < <?n. > " ; C 
HR* (3) 

1890 OUTPUT 10 : " 0 I ; F U " 

1908 FOP. 1 = 1 TO N 

1910 Z 1 1 I :• = I NT < < 0 < I ) -X 1> *83 > 

1 920 Z2 < I > = I NT < < P < I > - Y 1 ) *S4 > 

1930 OUTPUT 10 ; "PR" ,Z1 < I > , Z2<I > 

; "PD" 

1940 NEXT I 

1958 OUTPUT 10 ; "PIJ0, 908, 100, 980 

ll 



196 U 
1970 
1 980 

1 990 
2000 



OUTPUT 10 . "PIJ" 

OUTPUT 10 ,"81.22, 38" 

0I8P "INPUT THE LEGEND YOU 

*T NPMT IY E ’ 0 ‘ MHEM D0NE ” A label or legend can be entered here. 

If pL-0- them 2036 
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4. FIBER TEST MODULE (FTST) 



16 

20 

3 O 

4 0 
4 5 
50 



60 

70 

SO 



SO 
100 
1 1 0 
120 
1 30 
140 
150 



1 60 
179 

1 SO 
185 
1 90 
200 
210 



220 

230 

240 



250 

260 

270 

2 8 0 
290 

3O0 

310 



320 

330 

332 

350 

360 

370 

380 

3 9 0 

4 0 O 
410 



! F T S' T o 7 

< Jin nrgectte 

i 



JAN 



13 S 



The Fiber Test is known as "FTST". 



CLERP 



DISP " FT3T87 " 

DISP @ DISP "YOU MUST R'JU Cfl 
l.B BEFORE PLOT OR SAU-DAT 1 " 
DISP 0 DISP " USE THE FIJHCT I 
OH KEYS." 

SHOP T P < 450 > , L < 450 > , T 
SHORT D 1 '".450 :• . D2 < 450 > , G3< 450 
> • D < 4 5 O "> 

SHORT Z 1 <450 > , Z 2 4 5 O > 



"CALB" program must be run before the 
PLOT" or "SAV-DAT" routines. 



GOTO 120 
CLERP 

0 1 \ k E Y # 4 . " G R R F " G 0 S U B 2 7 0 0 

0 N K E Y # 1 , " I N P U T S " G 0 S I J B 2 3 0 

0 M K E Y # 2 , " P C Q U I P " G 0 S U B 6 O O 

OH KE't # 6.<"RD_DAT" GOSUB l 1 7 The function keys are defined. 

O 



OH KEY# 7 • " F LOT " GOSUB 1320 
OH KEY# 5, "SAV-DAT" GOSUB 20 



30 



OH KEY# 8 • "Cfl L B " GOSUB 1830 
OH KEY# 3, "AREA" GOSUB 2500 
KEY LABEL 
GOTO 120 

CLEAR A DISP "EHO FTST87" 0 
END 



! ft INPUTS **. 

CLEAR A i DISP" ENTER DATA R 
A T E < P T S MIN > " 



! I f) PUT R 
R=600 ! PRESET DATA 
DISP "ENTER ELAPSED 
OP 45 >" 

INPUT EO 

CLEAR ® DISP "ENTER 
m •••• fu l n ) " 

INPUT S 

CLEAR 0 DISP "ENTER 
GTH (*n'" 

INPUT G 1 
El=E0tl0OO 



RATE 

TIME C. 30 



SPEED Cm 



GAGE LEN 



The program loops here until a function is 
selected. 



The "INPUTS" routine is required first. 



The gage length is to be used in computing 
strain. 



CLEAR 

DISP 



DISP "ENTER ‘A 1 FROM LDCALB 
LISTING" 

INPUT C 3 

DISP "ENTER * B ’ FROM LDCALB" 
INPUT B 
CLEAR 70S 

OUTPUT 70S ; " A I 1 0 " 
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The excitation voltage is read and the 
calibration adjusted to compensate for any 
voltage drifL 



420 


ENTER 709 , 0 




430 


OISP "EXCITATION VOLTAGE IS" 
• 1 1 


4 40 


• v 

A=ABS< C3' > 




450 


A = I NT i At 1 000 > / 1 000 




460. 


DISP " A = " ; A 




4 70 


CUSP 




475 


OISP "SPEED = " ■ S : " N M 


/MIN" 


4 o 0 


COPY 




490 


CLEAR '2 OISP "ENTER 
CE COEFFEC I ENTS " 


COMPLIAN 


500 


OISP " E= " 




510 


INPUT E 




2 O 


OISP " F = " 




530 


INPUT F 




54 0 


DISP "6= " 




550 


INPUT G 




5 6 O 


PRINT " E = ", E G PRINT 


" F = " ; F 


570 


PRINT " G = " ; G 




5S0 


CLEAR 




530 


RETURN 




6 0 0 


i 




6 1 0 


! nr* hC QUIP i.r 




620 


IF E 0 > 0 THEN 660 




630 


BEEP (2 CLEAR 




640 


OISP "'INPUTS' MUST 
IRST i " 


BE PUN F 


650 


RETURN 




6 6 0 


CLEAR 




690 


BEEP 




700 


FOP 1 = 1 TO 7 12 OISP 


G NEXT I 


719 


OISP "CHECK IEEE BUTTON,. NOU 




NT SAMPLE AND BURN" 




720 


OISP "HIT CONT WHEN 


READY" 


722 


PAUSE 




723 


CLEAR 




724 


OISP "FINDING ZERO LOAD" 


725 


OUTPUT 704 " K 1 3 • " S 




730 


OUTPUT 709 ; " ARVR 1 VT3VD5VA0A 




E2 " 




731 


OUTPUT 709 : " A I O " 




732 


ENTER 70S ; L 




733 


P=A*L+B 




734 


IF P<0 THEN GOTO 755 




735 


OUTPUT 704 "K2" 




7 36 


WAIT 50O 




737 


OUTPUT 704 "KO" 




733 


GOTO 731 




750 


CLEAR G BEEP G OISP 
UP" 


"MOVING 


755 


OUTPUT 704 ; " K 3 " 




756 


OUTPUT 709 ; "RIO" 




757 


ENTER 709 ; L 




759 


F’ = R*L + B 





The compliance coefficients are entered so 
that a matching compliance curve may be 
generated and then subtracted from the total 
deformation, yielding the true displacement 
of the sample. 



The routine ends and returns to the menu. 

The "ACQUIR" routine is the actual data 
acquisition routine. 

The program checks to see that the inputs 
have been entered and provides prompts if no 
input data can be found. 



The OUTPUT ”704” codes provide for 
automatic control of the Instron tester. 
Initially the cross head is adjusted to provide 
for slight tension in the fiber, after which 
the timer will start and data collection will 
begin. 
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769 IF PC 06 THEN GOTO 75b 

76E OUTPUT 704 , "K21" The voltage from the load cell is converted 

765 BEEP 0 OISP "TIMER ON" to physical units of grams. Here the tension 

76P ON TIMER# l.-El GOTO 830 level is adjusted. 

770 1=1 

780 OUTPUT 70S 1 "RIO" 

790 ENTER 709 ; LCI> 

800 1=1+1 
810 WRIT 23 
P» 2 0 uOTO r' 8 0 
830 OFF TIMER# 1 
835 OUTPUT 704 ."KO" 

840 1=1-1 

8 50 WRIT 5 O 
860 CLERR 

8 70 FOR T=l TO 7 w OISP 8 NEXT 1 
880 BEEP (? CLERR l? OISP "CHECK E 
XTENSION HIT CONT" 

880 PR USE 



900 OUTPUT 704 ; " K 1 " 

910 N= I 

920 OISP ND'OATR POINTS" 

930 OISP "COMPUTING " 

940 FOR 1=1 TO N 
950 PCI ;=R*LT I >+B 
960 NEXT I 

9 70 P 2 = P C 1 > i? P1=PC1> 

980 FOR 1=2 TO N 

990 IF PC I >>P2 THEN K=I 

1 000 IF P < I > > P2 THEN P2 = P C I > 

101O IF PC I > CPI THEN P1=P CD- 

1030 NEXT 1 

1068 OISP " P < m ax > = " ; P2 After a quick analysis of some key data 

1070 D 4 = S t C K ✓ R > 4 f 1 G 1 > 1 1 00 points arc printed out to provide a indication 

1080 OISP " 0 X " D 4 > K of the nature of the data gathered. 

1090 PRINT " P C M R X > = " ; P2 
1095 PRINT "OX = ";D4,K 
1150 RETURN 



1 160 
1 170 
1180 
1 190 
1195 



1 200 
1210 
1 220 
1230 
1240 
1250 
1260 
1270 
1 280 
1290 
1 300 



! t* RERD ORTR ** 

CLERR 

OISP "ENTER FILE NRME" 

OISP "YOU MUST HAVE INPUTS 
BEFORE YOU CRN PLOT RERD DR 
Tfl" 

INPUT N$ 

RSSIGN# 1 TO NT 
RERD# 1 : M'E0iS.'Gl)K 

FOP 1=1 TO N 
RERD# 1 ; PCDDOCD 
NEXT I 

ASSIGN# l tq * 

Y2=P C K > R X 2 = D ( K ) 

OISP " ORTR HRS BEEN RERD" 

BEEP 

RETURN 



The sub routine returns to the main loop. 

The "RD-DATA” routine allows you to read 
in a data file that had been stored earlier and 
then plot results. 

It reads: Number of data points, N 

Elapsed time of test, EO 

Speed of test, S 

Gage length tested, G1 

Point of Max load, K 

Load and Strain, P(l), D(I) 
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I 

I 



** PLOTTER ** 



Program returns to the menu. 



1310 
1 320 



133 O 
1 340 


C L E R R 
CUSP " 

c «» 


ENTER 


Y NRX 


FOR Y HXI 


1350 
1 360 


INPUT 
DISP " 

it 


ENTER 


Y MIN 


FOR Y RXI 


1 3 7 O 
1330 


INPUT 
DISP " 

C M 


Y 1 

ENTER 


X MRX 


FOP X RXI 


1330 
1 4 0 0 


INPUT 
DISP " 

c }l 


V 7 

ENTER 


X MIN 


FOR X RXI 


1410 


INPUT 


XI 







The "Plotter" routine is the same as used in 
"MCTST". 



1420 S 3 = 1 0 0 0 /' < X - X 1 > 

1430 ' 34=1 0 0 0 ' < V 2 - Y 1 ) 

14ay PRINTER 13 10 
1450 CONTROL 105 • 43 

1 4b0 OUTPU T 10 ; " I N , SR 1 ; I P240O , 1 

60 O 3300 . 6300 ; " 

1470 OUTPUT 10 "300,1000/0,1000 



1 430 
1 4 90 



1 5 O 0 
1510 
1 520 
1530 
1540 



1 550 
1560 
1570 
1580 
1 590 
1 600 

1610 
1620 
1 630 

1646 



1 650 
I860 
1665 
1670 
1 680 
1690 



OUTPUT 10 . " Pi I© , 0 P 0 0 . 1 000, 1 
0 00, 1 0 0 0 , 1 0 0 0 , 0 . 0 , 0 p IJ " 
OUTPUT 10 "310 . 2, 0 3; TL1 . 5 

, 0" 

FOR 1=0 TO 5 

V = I NT < Y 1 + < Y 2 - Y 1 • ✓ 5*1) 

Y4 = TNT (. < Y-vi >*34> 

0 U T R IJ T 1 0 , "PH O , " , Y 4 , " Y T " 
OUTPUT 1 0 "CR-5, -O . 07.; LB" ; 

Y ; C H F: t < 3 ) 

NEXT I 

FOP 1=1 TO 5 
X = X 1 + ( X2-X 1 >.''5*1 
X 4 = I N T < < X - X 1 > 1 33) 

OUTPUT 10 " PR " ; X4 , " , O ; XT .; " 

OUTPUT 10 , "CP— 1 .3 , -1;LB" .; X 
; CHPf <3'.' 

NEXT I 

OUTPUT 10 "31.30, 42" 

OUTPU T 10 ? " PR220 • 0 ; CP 6 , -2 . 
3 .; LBSTRfi I N • V. > “ , C H P $ < 3 ) 
OUTPUT 10 , "PRO, 460.; OIO, 1 ; C 
P-2 6, 2. 60; LB FORCE < =rn ) " ; C 
H R f < 3 ) 

OUTPUT 10 "01; PU " 

FOR 1=1 TO N 
D- I ) =3*1 / R 

Z 1< I ) = I NT < <. 0 < I > -X 1 > * 33 > 

Z2 v I ) = I NT < (F'CI)-Yl ) *34 ) 

OUTPUT 10 .; "Pfl" , Z1 < I > , Z2CI ) 
. " PD '• 



Here, D(I) represents the strain while the 
load is at P(I). 



75 



ireo next i 

1 7 1 0 OUTPUT 10 . " PI 10 , 900, 1 00 , 900 

ii 



1720 OUTPUT 10 ,"plj" 

1730 OUTPUT 10 ■ "PI . 22 , . 38" 

1740 DISP 8 CUSP "ENTER THE LEGE 
NO YOU HR NT " 

1750 DISP "ENTER R 'O' WHEN DONE 

ll 



1760 INPUT Pi 



17 70 
1 780 



IF P i = " © " THEN 1800 
OUTPUT 10 ; " CP LB " Pi ; CHPi < 



Label the graph for identification. 



1 790 
1 800 
1818 
1 820 
183 0 
1 840 
1 850 
1 86u 
1 870 

1880 
1 890 
l 9 0 0 
1910 
1920 



193 0 
1 940 
1 950 
1 960 
1970 
1 980 



1 990 
2000 
2010 
2815 



GOTO 1^40 
PRINTER IS 2 
RETURN 

i 

' M * CfiLB ttt 
IF F < > 0 THEN GOTO 1870 
CLERP i? DISP "HEED INPUTS'" 
RETURN 

CLEAR 0 BEEP 6 DISP " DRT A C 
OH VERS ION TO PHYSICRL UNIT" 
FOR 1=1 TO K 
Did >=STI/P 

D2< I > = G 4 P < I )^2 + F*P(I >+E 
D3< I > =D 1 < I >-D2< I > 

DC I >=D3< I ')'■ Gli 100 ! X STRRI 
N 

NEXT I 

FOR I = K + 1 TO N 
D1C I > = S 4 l F 
D2< I > = D 2 C K > 

D3C I >=01 C n-C'2< I > 

D < I > = D 3 < I > / G 1 i 1 0 0 ! STRRI 
N 

NEXT I 

PRINT " P C M R X > = " : POO ; K 

PRINT "D*: =" iDCfO 

CLERR R DISP "CRLB COMPLETE 



The routine returns to the menu. 

"CALB” is the calibration routine where die 
final output is true displacement and strain 
of the fiber sample. 

D1 is the total displacement. D2 is the 
displacement of the system as a function of 
the load, PI, calculated using the coefficients 
E, F, and G. 

D3 = D1 -D2, the true displacement of the 
fiber sample. 

D is the strain in percent . 



202O 

2030 

204O 

2050 



206 O 
20 70 

2080 
2090 
2 1 GO 
2110 
2120 
2130 
2140 
2150 



RETURN 

! *** STORE ADJUSTED DRTR t 
4: 

IF DC 30>O0 THEN 20 7 O 

BEEP 0 DISP " NO DRTR PRESE 

NT" 

RETURN 

CLERP l? BEEP 12 DISP "INSERT 
DRTR CARTRIDGE" 

DISP "ENTER FILE NRME . " 
INPUT Hi 

CLERR I? DISP "STORING DRTR" 

N = 2 4 CN+l >4:8+50 

H=INT(P1/'256> + 1 

CRERTE N $ , H • 2 5 6 

ASSIGN* 1 TO Hi 

PRINT# 1 .; N , E 0 • S .• G 1 , K 



Return to the menu. 



This sub routine is used to store die adjusted 
data calculated in the CALB sub routine. 
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21-50 
2170 
2130 
2130 
22O0 
2210 
2220 
250O 
25 1 O 
2520 



2 5 3 0 
2540 
2550 
25-50 
2570 
2530 
2530 

2535 

2600 



2610 



2700 



FOR 1=1 T0 M 
PRINT# 1 } P<I),D<I> 

NEXT I 

ftSSIGN# 1 TO % 

CLEftR i? DIOR " DftTft STORED" 

BEEP 

RETURN 



' H ftp Eft ** 

CLEftR 0 D I S P " COMPU T I NG HR E Return to menu 
ft UNDER P/D CURVE " 

T=U AREA calculates the are under the load 

FOR 1 = 1 T 0 K - 1 versus displacement curve, 

ft 1 = •' P ( I > +P< I + 1 > > /2 
ft 2 = ft 1 1 ( D < I + 1 > - D I > ) 

T = T + ft 2 

NEXT I 

PRINT "ftREft UNDER THE P/0 C 
URWE=" ; T; " qmJj" 

PRINT 

IF N T < > " " THEN PRINT ; " 



" • NT 

CLEftR ft DISP "ftREft COMPUTED 
" ft RETURN 



i 



2710 ! * ID ISP CRT GRAPH*-* 

2720 PRINT " " NT ft 

PRINT 

2 7 3 0 C L E ft R ft D I S P “ I N P U T P < M ft X > 
FOR V ft X IS" 

2740 INPUT P2 
2750 P1=0 
2 7 6 0 G C L E ft R 
2730 GRftPH 
2 730 L Cl I P O 

2800 SCHLE -48, 208.. -36/ 156 
281 U X ft X I S O .. 32 , 0 , 1 60 
2820 YftXIS 0.. 28.. 0 .140 
2830 FEN 1 
2840 PE HUP 

2 8 50 81 = 1 6 0 < E 0 6 0 > / S 

2860 82= 1 40/P2 

2870 FOR T= 1 TO 5 

2880 X = S * T / 1 0 

2830 X= I NT (X* 1 000) /I 000 

2300 X2=X*S1 

2910 HONE X 2 , -10 

2920 LftBEL VflL*<X> 

2930 NEXT T 
2940 FOR T=0 TO 5 
2950 p = P 2 / 5 * T 
2960 P3=P*S2 
2970 MODE -23- P3 
2980 LftBEL VftLTCP) 

2930 NEXT T 
3090 MOVE 40,-23 
3010 LftBEL "DISPCMM) " 



This sub routine provides for a graphical 
output on the CRT screen 



77 



3020 
3 0 3 O 
304 0 
3050 
3060 
3070 
3080 
3090 
3 1 00 
3110 
3 1 20 
3 1 30 
3 1 40 
3 1 50 
3 2 0 0 
3210 
3220 
3230 



3240 



3250 

3251 
.7252 

3253 

3254 
3260 
3270 
3280 
3 2 9 0 
3300 
33 1 0 
3320 
333 0 
3340 
3350 
3360 



3370 



DEG 

MOVE -29 , 25 
LDIP 90 

LftBEL "FORCE (?n.) " 

LDIR 0 

FOR 1=1 TO M ' 

D< I >=S*< I / R 
Z 1 < I > = D < I > Z S 1 
Z2< I >=P< I >*82 
PLOT 21 < I > , Z2< I > 

NEXT I 
COPY 
CLEAR 
RETURN 

! tZO I AMETER M 

! DIAMETER FROGRAM 
CLEAR 

i U=HAVE LENGTH, N T = F I B E R 
ESI GNAT I ON, Cl-OELTAC AND 
3=SEM EQUIVELEHT VALUE 
DISP "THIS PROGRAM FINDS T 
E FINAL DIAMETER" 

W= . 0000006328 
N= 1 
C=1 15 
L = 6 9 8 1 

DISP "ENTER FIBER DESIGN" 
INPUT NT 

PRINT "FIBER NO ";N$ 

FOP 1=1 TO N 

DISP "ENTER MICRO READING" 
INPUT Cl 
Q = U 4 2 1 L ■ ' C 1 + C > 

0= I NT f. 0* 1 E 1 3 > ••••• 10000000 
PRINT "NPSLT VALUE = " G 
! CONVERT TO SEM VALUE 
Ql = . 8 9 1 Q + .3614 
0 1 = I NT < Q 1 * 1 0000000 ) / 1 000000 
0 

PRINT "FIBER DIAMETER- ";Q1 
}” MICRO 

NS" 



This sub routine is used when determining 
fiber diameter, it provides for the input of 
the fiber designation and micrometer 
readings. 

The fiber diameter is determined based on an 
average of the micrometer readings and 
corrected in accordance with the SEM cal 
curve. 
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5. LOAD CELL TEST MODULE (LDTST) 



i o 


it LDTST M- 




20 


' TAKES 25 READINGS 


FROM 




fiDCELL AMD PRINTS ftV 


ERAGE 


30 


C L E A P 




4 0 


CUSP •‘EATER A " 




J0 


IMF UT Hi 




Mj 


C'ISP "ENTER E" 




70 


INPUT B 1 




8 0 


CLEAR 




30 


CLEAR 702 




1 0 0 


OUTPUT 70? • " A I 19 " 




1 10 


ENTER 70? ; V 




120 


A = A 1 0 




1 25 


B = B1 




1 30 


P 1 = 0 




1 4 0 


FOR 1=1 TO 25 




1 50 


OUTPUT 70S ; "RIO" 




1 60 


ENTER 70? , L 




1 70 


P = L* A + E. 




1 75 


P= I NT T P T 1 0O0 > "" 1 060 




1 SO 


OISP F‘, 




1 ? 0 


P1=P1+P 




2 0 0 


ne::t i 




210 


PRINT “VC out > = " , v' 




220 


PRINT "LOAD = M iPl/25 


.• " 9 In 


2 3 O 


BEEP 




240 


PRINT 




2 50 


OISP "SET LOAD. HIT 


CO NT" 


2 b 0 


PAUSE 




2 7 O 


PRINT 




2S0 


G 0 T 0 :? 0 





LO "LDTST" is a short utility to test ihc load 
cell for proper operation and reproducibility. 



Load cell output voltage is recorded 25 times 
then converted to grams. 

The average load is printed. 



The program pauses and a load level may be 
changed. 



Press continue after installing new load. 



Press "pause” or "reset" to leave program 
loop. 
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6. LOAD STORAGE MODULE (LOADSTOR) 



13 ' ** LOST OF’ U 

1 LOROS ROD STORES OhTh 
TAPES CREATED WITH FTS 
MCTST 

« JIM N A G E 0 T T E OCT 27 
DIM P 1C 450 > .• DM50 > 

G 0 S U B 5 0 6 0 
CLEAR 

"PD-C" GOSUB 
" SA'J-C " GOSUB 
"FD-F" GOSUB 
"SHV.'-F" GOSUE' 
" EHD" GOTO 15 



30 
40 
4 5 
50 
so 
70 
80 
90 
1 OO 
1 10 
1 20 
150 
1 OOO 
1010 
1 020 

1 030 
1 040 
1050 
1 060 
1 070 
1 08O 
1 090 
1 100 
1110 
1 120 
1 1 30 
1140 
1 150 
1160 

1 170 
2O0O 
2010 

2020 

2 0 3 0 
2040 
2050 
206O 
2O70 
2030 
2890 
2100 
2110 
2120 
2130 
2140 
2 1 50 
2160 
2170 
2130 



OH 

OH 

OH 

OH 

OH 

KEY 



1 



KE'r # 

KEY# 
KEY# 
KEY# 
KEY# 8 
LABEL 
68 



4. 



FILE 
T AND 
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1 0 0 0 
2 O O O 
3 0 0 0 
4 000 
O 



GOTO 
EHD 

! *:* PD-C ** 

CLEAR i? BEEP 

CUSP u READ COMPLIANCE 



FILE 



DI3P 
OISP 
DISP 
IHPUT 
CLEAR 
ASSIGN# 

READ# 1 
FOR 1=1 
READ# 1 
WENT I 
ASSIGN# 

CLEAR 

CUSP "DATA 
EEEP 
RETURN 

• ** SAVE C ** 
CLEAR A DISP '* 

CE FILE" 

DISP 

DISP "INSERT DATA 
BEEP 
DISP 
INPUT 
CLEAR 

DISP "SOWING DATA" 
M = < N + 1 >*8 + 50 
H= I NT < M ' 256 > + ! 
CREATE N$ > H / 256 



INSEFT DATA TAPE" 
ENTER FILE NAME" 

NT 

S DISP "READING DATA' 
1 TO NT 
; H . EO / S 
TO N 
s P 1 < I > 

1 TO * 

READ" 



AVE COMPLIAN 



"LOADSTOR" is a utility program to allow 
for quick storage or access to data files. This 
is useful for the transfer of files. 



"RD-C" and "SAV-C” read and save the 
"MCTST" files. 



"RD-F" and "SAV-F" read and save the 
"FTST" files. 



After loading a data file, you may press 
"END" and then examine or list the data in 
the command mode. Refer to the print out 
for variable names. When action complete, 
re-enter the program without losing the data 
by entering "cont 50". 



ENTER 

NT 



TAPE' 



FILE NAME" 



ASSIGN# 

PRINT# 

FOR 1=1 

PRINT# 

NEXT I 

ASSIGN# 

CLEAR 



1 TO 
; N 
TO N 
; PI < I > 



NT 

E0, S 



1 TO * 
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DISP "DATA SAWED 
BEEF' 

RETURN 



** F'D-F %% 



7. @ 1 @ CLEAR 6 BEE P 

3020 DISP “ READ FIBER TEST FILE 



3 0 3 0 D I S F 

3040 DISP "INSERT DATA TAPE” 

3050 DISP "EHTEF FILE NAME" 

3060 INPUT Nf 

3O70 CLEAR 6 DISP "READING DATA" 

3080 ASSIGN# 1 Tfi HT- 

3090 READ# 1 ; N,E0,S,G1.K 

31O0 FOR 1=1 TO U 

3110 READ# 1 ; PI ( I > , D< I ) 

3120 NEXT I 

3130 ASSIGN# 1 TO % 

3140 CLEAR 

3150 DISP "DATA READ" 

3160 BEEP 
3170 RETURN 
40O0 ! ** SAVE F t t 
4010 CLEAR R DISP "SAVE FIBER TE 
ST FILE" 

4020 DISP 

4030 DISP "INSERT DATA TAPE" 

4040 BEEP 

4050 DISP "ENTER FILE NAME" 

4O60 INPUT Nf 
4O70 C L E A R 

4030 DISP "SAWING DATA" 

4030 N=2*<N+i:-*S + 50 
4100 H=INT' fV256> + l 
4110 CREATE N $ .• H ■ 256 
4120 ASSIGN# 1 TO N* 

4130 PRINT# 1 N. EO , S . G 1 , K 

4140 FOR 1=1 TO N 

4150 PRINT# 1 ; F'1(I),0(I) 

4160 NEXT I 

4170 ASSIGN# 1 TO * 

4180 CLEAR 

4130 DISP "DATA SAWED" 

42O0 BEEP 

4210 RETURN 

5000 ! ** INFO * * 

5010 PRINT "AFTER 'PAUSE' OP ' EN 
D' YOU MAY CONTINUE WITHOUT 
DATA LOSS BY TYPING 
5020 PRINT ’"CONT 50’" 

5030 PRINT 

5840 PRINT "C FILES CONTAIN N,EO 
, S AND P1<I>" 

505O PRINT 

5060 PRINT "F FILES CONTAIN N,E0 
, S , G 1 , K AND P 1 ( I > / D < I > " 

507O PRINT 8 PRINT 



II 
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APPENDIX C. DETERMINING BUNDLE PRELOAD 



Before composite preload level can be determined a mathematical 
model, based on individual fiber data, must be developed to study the effects 
of gage length on the strength distribution. Previous studies have provided 
conclusive results that composite strength is a function of gage length. In 
order to predict the behavior of a 254 millimeter strand, the results of the 
fiber tests, (gage length = 50 millimeter) must be examined and equated to a 
sample that is 254 millimeters long. 

Results from a previous study [ Ref. 7], which examined the affects of 
proof testing Kevlar-49 at elevated temperatures will be used to verify the 
mathematical model and to establish an estimated preload level. 



1. MATHEMATICAL CORRECTION FOR GAGE LENGTH 



FI = 1 - exp{-(a/(3i) a i} for gage length = \\ (4) 



R1 = exp {- a / P i ) a 



( 5 ) 



l 2 /ll = n 

Rt = n Ri R 2 — • = (Ri) n 



for gage length = I 2 



= [exp {-o/pi) a i}} n 



= exp { -n ( a / Pi ) a i } = exp { - n(a /P2 ) “2 } 

= exp { -n ( o a i / Pi a i ) } = exp { -n(c a 2 / p2 a 2)} 



a a l - a a2 . . ai = a2 
n / Pi«! = 1 / p 2 «2 

( Pi / P 2 ) a i = n = h /l 2 
( pi / p 2 ) = l 2 ( 1/a l) / lid/ai) 
(p2/pl)= Oi/l 2 )(!/«!) 



Results: 



P2 = Pi ( *i / h) 



( 6 ) 
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Mathematically, equation (6) provides for the beta adjustment for gage 
length . 

2. KEVLAR -49 MODEL 

Based on original test data. Table VII., the beta for strain (3e can be 
corrected for various gage lengths. 



TABLE VII. ORIGINAL TEST DATA 



Material 


G.L.(mm) 


a 


P(gm) 


(gm/8) 


Kevlar -49 


50 


6.7 


38.9 


1.3x103 


Kevlar (25 °q 


254 


25.9 


9.74x103 


4.0x105 


Kevlar (700Q 


254 




9.50x103 


3.0x105 



Using the test data for a 50 mm. fiber and substituting (3 e for (3j in 
equation (6), where; 

Pe = Jl 
E 

Pe2 = Pel Ol / 1 2 ) 

Correction to gage lengths can be made, see Table VIII. 

TABLE Vm. CORRECTED pe OF KEVLAR - 49 FOR 
VARIOUS GAGE LENGTHS 5. 



Delta mm 
(ineffective length) 


5 


.1 .5 1 


10 


50 


254 


Beta of Strain 


Pe 


.076 .059 .054 


.038 


.030 


.024 



Using the reults from a proof test of a Kevlar-49 strand (gage length = 
254 mm) at 70°C the number of individual fiber breaks can be determined. 
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pe = p L /£ 



Pe = 9.50xl0 3 gms / 3.0xl0 5 gms/e 
pe= 3.1x10-2 

Note : This result is somewere between 5=10 mm. and 50 mm. (Table 
VIII.) and correcting for gage length , equation (6), 5 = 50 mm. . 

The percent of fiber failures F as a function of a possible ineffective gage 
length is determined by equation (4). 



F = 1 -exp {-( P £GL /Pe70c ) a l} 



TABLE IX. PERCENT FIBER BREAKS AS A FUNCTION OF 
OF INEFFECTIVE GAGE LENGTH 70°C 



5 


Pe 


F 


#of fiber breaks 
(267 per bundle) 


.1 


.076 


.0001 


.03 


.5 


.059 


.005 


1 


1.0 


.054 


.01 


3 


10. 


.038 


.1 


27 


50 


.030 


.39 


104 


100 


.027 


.63 


168 


254 


.024 


.99 


267 



Based on the fact that the strand was heated to 70°C it would not be 
unreasonable to expect that the ineffective length 5 approaches 50 mm. 
Therefore from Table IX, 39 percent of the fibers broke. 
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3. GRAPHITE TESTING 

The model was then examined for graphite fibers based on known test 
data from individual fiber testing conducted during this study. For gage 
length = 50 mm; 

Pl = 15.5 gmf a = 4.0 p e = 1.72 x lO'^mm/mm E= 900 gm/e 

Using equation (6), Pl = 13.0gmf, for gage length = 100 mm. Now 
applying this information to a dry bundle consisting of 3000 fibers 254 mm. 
long the number of fiber failures can be predicted as a function of 
displacement using equation (4), (Table X.) 

F = 1 - exp{ - (£/pe) a e } 

P e 50 mm = -0145 mm/mm 



TABLE X. PERCENT FIBER BREAKS AS A FUNCTION OF 
STRAIN FOR GAGE LENGTH OF 254 MM. 



£ 


AL (mm) 


F 


No. of Fil 


.001 


.254 


.00006 


<1 


.002 


.508 


.0009 


3 


.004 


1.016 


.015 


44 


.006 


1.524 


.07 


214 


.008 


2.032 


.21 


626 


.010 


2.54 


.46 


1306 


.012 


3.04 


.69 


2083 


.014 


3.56 


.89 


2666 


.016 


4.06 


.98 


2929 


.018 


4.57 


.99 


2992 


.020 


5.08 


.999 


2999 



For example, based on what is know of the successful Kevlar test at 70 
°C (8=50 mm) it is desired to break 40 percent of the fibers during 
preloading. So now, the graphite bundle would be loaded until delta length 
equals 2.5 mm. (Table X). 
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Load KG 



These results can be graphical confirmed by reviewing a representative 
load-displacement plot (Figure 18) of a dry bundle (gage length = 254 mm.). 




Figure 18. Load-Displacement Plot 
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APPENDIX D. TABULATED TEST RESULTS 



TABLE XI. SINGLE GRAPHITE FIBER 

DIAMETER RESULTS (MICRONS) 

008 SERIES 019 SERIES 



6.707 


7.282 


6.208 


7.291 


6.784 


7.283 


6.435 


7.300 


6.797 


7.285 


6.515 


7.310 


6.800 


7.285 


6.593 


7.314 


6.815 


7.290 


6.730 


7.316 


6.852 


7.294 


6.798 


7.333 


6.866 


7.306 


6.804 


7.335 


6.898 


7.307 


6.815 


7.338 


6.918 


7.318 


6.825 


7.345 


6.940 


7.321 


6.889 


7.350 


6.964 


7.322 


6.915 


7.350 


6.965 


7.324 


6.945 


7.351 


6.973 


7.330 


6.962 


7.364 


6.983 


7.335 


6.996 


7.372 


6.997 


7.336 


6.998 


7.372 


6.997 


7.337 


7.015 


7.380 


7.020 


7.341 


7.030 


7.389 


7.039 


7.358 


7.042 


7.395 


7.045 


7.359 


7.054 


7.426 


7.060 


7.370 


7.062 


7.430 


7.066 


7.376 


7.062 


7.440 


7.084 


7.377 


7.070 


7.447 


7.104 


7.377 


7.101 


7.453 


7.115 


7.378 


7.137 


7.454 


7.140 


7.388 


7.145 


7.458 


7.141 


7.397 


7.147 


7.480 


7.142 


7.403 


7.178 


7.482 


7.179 


7.420 


7.180 


7.507 


7.184 


7.425 


7.197 


7.516 


7.209 


7.445 


7.199 


7.551 


7.227 


7.456 


7.204 


7.569 
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7.236 


7.473 


7.206 


7.585 


7.241 


7.573 


7.208 


7.586 


7.243 


7.576 


7.212 


7.589 


7.255 


7.650 


7.212 


7.600 


7.258 


7.682 


7.216 


7.623 


7.270 


7.760 


7.219 


7.670 


7.275 


7.836 


7.223 


7.738 






7.227 


7.792 






7.243 


7.825 






7.289 


7.871 
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TABLE XII. SINGLE GRAPHITE FIBER 
FAILURE LOAD RESULTS 



019-N 


019-1 


008-N 


008-1 


8.096 


4.303 


10.592 


8.022 


9.284 


5.015 


10.654 


9.255 


9.734 


6.819 


10.776 


9.653 


10.219 


7.368 


11.451 


10.707 


10.342 


8.111 


11.483 


11.018 


10.772 


9.899 


12.043 


11.157 


11.686 


10.324 


12.123 


1 1.243 


12.505 


10.636 


12.141 


11.622 


12.690 


11.751 


12.252 


1 1.781 


12.768 


11.770 


12.398 


1 1.804 


13.180 


12.440 


12.439 


12.499 


13.196 


12.448 


12.664 


12.623 


13.457 


12.667 


12.858 


13.071 


13.762 


12.904 


13.604 


13.157 


13.872 


13.703 


14.440 


13.297 


13.980 


13.705 


14.614 


13.630 


14.271 


14.265 


15.107 


14.140 


14.528 


14.502 


15.209 


14.812 


14.678 


14.528 


15.241 


14.914 


14.767 


14.811 


15.372 


14.936 


15.033 


15.565 


15.567 


16.1 16 


15.521 


16.957 


15.851 


16.369 


15.522 


17.073 


16.270 


16.633 


15.615 


17.726 


16.272 


16.941 


15.956 


17.804 


16.664 


17.560 


16.145 


18.259 


17.352 


17.607 


16.242 


18.801 


17.360 


17.923 


16.369 


18.956 


18.009 


18.992 


17.545 


19.392 


18.064 


20.603 


18.081 


19.556 


18.586 


20.91 1 


19.962 


19.610 


19.544 


21.572 


20.352 


19.915 


20.490 


25.442 


20.627 


N = NPSIT Results 


I = INSTRON Results 
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1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 



TABLE XIII. FAILURE LOAD RESULTS OF 
SERIES 019 GRAPHITE 



LOAD 


DIAMETER 


TRIAL 


LOAD 


DIAMETER 


(gms) 


(microns) 


(N) 


(gms) 


(microns) 


3.447 


7.219 


1 


7.284 


7.585 


4.303 


7.440 


2 


8.096 


6.815 


5.015 


7.426 


3 


9.284 


6.800 


6.819 


7.430 


4 


9.734 


7.242 


7.368 


6.992 


5 


10.219 


7.350 


8.111 


6.998 


6 


10.342 


6.515 


9.899 


7.364 


7 


10.772 


7.344 


10.324 


7.145 


8 


11.686 


6.593 


10.636 


7.054 


9 


12.505 


7.227 


11.751 


7.335 


10 


12.690 


7.482 


11.770 


7.333 


11 


12.768 


7.516 


12.440 


6.730 


12 


13.180 


7.589 


12.448 


7.030 


13 


13.196 


7.458 


12.667 


7.042 


14 


13.457 


7.453 


12.904 


7.147 


15 


13.762 


7.507 


13.703 


6.804 


16 


13.872 


7.600 


13.705 


7.289 


17 


13.980 


7.447 


14.265 


7.197 


18 


14.271 


7.015 


14.502 


6.825 


19 


14.528 


7.454 


14.528 


7.825 


20 


14.678 


7.380 


14.811 


7.569 


21 


14.767 


7.206 


15.565 


6.945 


22 


15.003 


7.551 


16.957 


7.062 


23 


15.521 


7.137 


17.073 


7.310 


24 


15.522 


6.889 


17.726 


7.291 


25 


15.615 


7.350 


17.804 


7.316 


26 


15.596 


7.623 


18.259 


7.738 


27 


16.145 


7.101 


18.801 


7.351 


28 


16.242 


7.208 


18.956 


7.372 


29 


16.369 


7.872 


19.392 


7.199 


30 


17.545 


7.212 


19.556 


7.216 


31 


18.081 


7.338 


19.610 


7.070 


32 


19.962 


7.389 


19.915 


6.962 


33 


20.352 


7.586 






34 


20.627 


7.792 
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1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 



TABLE 



LOAD 

(gms) 

5.504 

8.022 

9.255 

9.356 

9.653 

10.707 

11.018 

11.157 
11.243 
11.622 
11.781 
11.804 
12.499 
12.623 
13.071 

13.157 
13.297 
13.630 
14.140 
14.812 
14.914 
14.963 
16.116 
16.369 
16.633 
16.941 
17.560 
17.607 
17.923 
18.992 
20.603 
20.911 
21.572 
25.442 



XIV. FAILURE LOAD RESULTS OF 
SERIES 008 GRAPHITE 



DIAMETER TRIAL 



(microns) (N) 

7.283 1 

7.388 2 

6.852 3 

6.815 4 

6.965 5 

7.241 6 

6.918 7 

7.066 8 

7.285 9 

7.084 10 

7.060 11 

7.330 12 

7.377 13 

7.576 14 

7.358 15 

7.306 16 

7.285 17 

7.258 18 

7.224 19 

7.184 20 

7.045 21 

7.142 22 

6.973 23 

7.020 24 

6.944 25 

7.456 26 

7.397 27 

7.473 28 

7.335 29 

7.420 30 

6.997 31 

7.682 32 

7.642 33 

7.760 34 



LOAD DIAMETER 
(gms) (microns) 

10.234 6.832 

10.592 7.336 

10.654 6.707 

10.776 7.403 

11.451 6.964 

11.483 7.445 

12.043 7.236 

12.123 7.140 

12.141 7.039 

12.252 7.321 

12.398 7.282 

12.439 7.318 

12.664 7.141 

12.858 7.307 

13.604 7.337 

14.440 7.322 

14.614 7.275 

15.107 6.784 

15.209 7.294 

15.241 7.378 

15.372 7.115 

15.567 7.376 

15.851 7.270 

16.270 7.425 

16.272 7.179 

16.664 7.377 

17.352 7.359 

17.360 7.209 

17.360 7.243 

18.009 7.836 

18.064 6.997 

18.586 7.324 

19.544 6.800 

20.490 6.983 
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TABLE XV. COMPOSITE STRAND FAILURE 
LOAD RESULTS 



PRELOADED 

(KG) 



NON-PRELOADED * 
(KG) 



43.17 


46.76 


40.80 


47.29 


43.38 


46.79 


42.18 


47.54 


43.41 


46.95 


42.38 


47.69 


43.81 


47.28 


43.37 


48.05 


44.40 


47.46 


44.48 


48.31 


44.75 


47.84 


44.86 


48.44 


44.94 


47.97 


44.89 


49.08 


45.69 


49.02 


45.18 


50.35 


45.72 


49.91 


45.61 


50.60 


45.74 


49.99 


46.53 


50.80 


46.04 


50.17 


46.79 


51.91 


46.42 




47.15 


52.22 



* Preliminary data on AS4-Graphite Strand Intrinsic Strength, 

By: E.M. Wu, N.Q. Nguyen and G.W. Nypiuk, Lawrence Livermore 
National Labratory, Livermore, California, April, 1972 
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